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The gauge boson pair production processes 14^7, WW, WZ, and Z7 were studied using pp 
collisions corresponding to an integrated luminosity of approximately 14 pb"'^ at a center-of-mass 
energy of ^ = 1.8 TeV. Analysis of W"f production with subsequent W boson decay to iu {£ = e, /i) 
is reported, including a fit to the pr spectrum of the photons which loads to limits on anomalous 
WW^ couplings. A search for WW production with subsequent decay to livP [t = e, /x) is presented, 
leading to an upper limit on the WW production cross section and limits on anomalous WW^ and 
WWZ couplings. A search for high pT W bosons in WW and WZ production is described, where 
one W boson decays to an electron and a neutrino and the second W boson or the Z boson decays to 
two jets. A maximum likelihood fit to the pr spectrum of W bosons resulted in limits on anomalous 
WW-y and WWZ couplings. A combined fit to the three data sets which provided the tightest 
limits on anomalous WWy and WWZ couplings is also described. Limits on anomalous ZZ7 and 
Z77 couplings are presented from an analysis of the photon Et spectrum in Zf events in the decay 
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channels (ee, fijj,, and w) of the Z boson. 



I. INTRODUCTION 

Interactions between gauge bosons, the W boson, Z boson and photon, are a consequence of the non-AbeUan 
gauge symmetry of the Standard Model (SM). The gauge boson self- interactions are described by the trilinear gauge 
boson vertices and contribute to gauge boson pair production in pp collisions. The cross sections of these processes 
are relatively small within the SM. The inclusion of non-SM (anomalous) couplings at the trilinear gauge boson 
vertices enhances the production cross sections of gauge boson pairs, especially at large values of the gauge boson 
transverse momentum px, and at large values of the invariant mass of the gauge boson pair system. Observation of 
anomalous gauge boson pair production would indicate physics beyond the SM. Feynman diagrams for gauge boson 
pair production are shown in Figs. |^(a)-(c), where Vb,l^i and V2 are the W boson, the Z boson, or the photon. 
Figures |l|(a) and (b) are described by well-known couplings between the gauge bosons and quarks. Figure [l](c) shows 
the trilinear coupling diagram. Numerous phenomenological studies [Q of the characteristics of gauge boson self- 
interactions have been performed in anticipation of hadron and e^e~ collider experiments where direct measurements 
of the coupling parameters are possible by studying gauge boson pair production processes. 

This paper describes studies of gauge boson pair production and the corresponding trilinear gauge boson coupling 
parameters using data from y/s = 1.8 TeV pp collisions taken with the D0 detector during the 1992-1993 Tevatron 
collider run at the Fermi National Accelerator Laboratory. Four processes were studied: W"f production, where the 
W boson decayed into ev or fiv; W boson pair production, where both of the W bosons decayed into eiy or /ii/; 
WW and WZ production, where one W boson decayed into ei/ and the second boson decayed hadronically; and Zj 
production, where the Z boson decayed into e'^e^ , fi'^ fi^ , or i/v. 

This paper presents the details of analyses whose results have already been published (2|-|^. In addition, it presents 
limits on anomalous trilinear couplings from the combined W^, WW , and WZ analyses. 

A. WW^ and WW Z Trilinear Gauge Boson Couplings 

A formalism has been developed to describe the WWj and WWZ interactions for models beyond the SM using 
an effective Lagrangian. The WW-f and WWZ vertices that satisfy Lorentz invariance and conservation of C and 
P can be described by a Lagrangian with two overall coupling strengths gwwj = —e and gwwz = —ecotOw and 
six coupling parameters gY , ny and Ay, where V = 'y, Z . In the SM, AgY — gX — 1 = 0, Akv = Ky — 1 = 0, and 
\v = 0. The anomalous couplings are parameterized as form factors with a scale. A, in order to avoid unitarity 
violation of the gauge boson pair cross section at asymptotically high energies: e.g., A^(s) ~ A-y/(l + s/A^)^. The 
WWj and WWZ coupling parameters, in the static limit, are related to the magnetic dipole moments {i^w) £^nd 
electric quadrupole moments (Q^) of the W boson: = 2j^|-(l + K + X) and = -jp-{K - A) where e 

and Mw are the charge and the mass of the W boson. A more detailed discussion of the effective Lagrangian for the 
WW J and WWZ interactions can be found in Appendix 1. 

The Wj production process has the highest cross section among the gauge boson pair production processes at the 
Tevatron. Feynman diagrams for the qq' — > W^ process are obtained by substituting Vq = Vi = W and V2 — j in 
Figs. |l|(a)-(c). A delicate cancellation takes place between the amplitudes Q that correspond to the u and t-channel 
quark exchange, shown in Figs. |^(a) and (b), and s-channel production with a W boson as the mediating particle, 
shown in Fig. ^c). A W boson is identified by its leptonic decay products: a high px charged lepton £ {£ — e,/i); 
and large missing transverse energy (^t) due to the undetected neutrino. Single W boson production, followed by 
radiation of a photon from the charged lepton from the W boson decay, also contributes to the tiyj final state. This 
process is shown in Fig. ^ The photon from the radiative decay is preferentially emitted along the direction of the 
charged lepton; the process can be suppressed by imposing a minimum separation requirement, ARfy, between the 
charged lepton and the photon where ARg^ is the distance in pseudorapidity and azimuth. For ^/s = 1.8 TeV pp 
collisions, the predicted cross section times branching fraction for W ^ ev 01 jiv for W^ final states with photon 
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transverse energy > 10 GeV and ARi^ > 0.7 is 12.5 pb. Figure || shows the Et spectrum of photons for 

the SM and non-SM production processes predicted by the leading order theory ||l^]. The distributions for non-SM 
WW^ coupHngs exhibit a large increase in the cross section at high E^. The 14^7 process is sensitive only to WW^ 
couplings, not to WWZ couplings. It is more sensitive to than to Ak^, since the amplitudes related to and 
Ak^ are proportional to s and -s/l, respectively. 

Limits on the WWj trilinear couplings from studies of Wj production have been reported |ll|,|l^,^ previously 
by the UA2, CDF, and D0 Collaborations. In this paper, the results from the D0 Collaboration, with the most 
restrictive WWj limits, are presented in more detail than in the recent publication 

The process pp WW + X is predicted to have a cross section of 9.5 pb ||l^] at next-to-leading order at y/s — 1.8 
TeV. The Feynman diagrams for the WW production processes are obtained by substituting Vb = 7 or Z, Vi = W'^ 
and V2 = W~ in Figs. |l](a)-(c). Destructive interference, similar to that occuring in W^ production, occurs between 
the u and t-channel amplitudes and the s-channel amplitude with a photon ot a, Z boson as the mediating particle. 
The former processes contain the well-known couplings between the W bosons and quarks and the latter the WW^ 
and WWZ trilinear couplings. W boson pair production is sensitive to both of the WW^ and WWZ couplings. It is 
approximately a factor of two more sensitive to the WWZ couplings, due to the higher value of the overall coupling 
gwwz = —ecotOw, than to the WWj couplings with gww-r = — e and is therefore complementary to W^ production. 
The predicted cross section for WW production, as a function of anomalous coupling parameters A = A^ = A^ 
and An = Ak^ = Akz, where the WWj and WWZ couplings are assumed to be equal and A — 1000 GeV, is shown 
in Fig. |. 

The details of the recently published analysis |^ , in which an upper limit on the WW cross section was obtained 
from the observed number of signal events in the dilepton decay modes, are presented. The limit on the cross section 
was translated into the limits on anomalous coupling parameters. 

For pp collisions at = 1.8 TeV, the cross section predicted |l^ at next-to- leading order for WZ production 
is 2.5 pb. The Feynman diagrams for WZ production are obtained by substituting Vq = Vi = W and V2 = Z in 
Figs. |l(a)-(c). While WW production is sensitive to both the WW^ and WWZ coupHngs, WZ production depends 
only on the WWZ couplings. The WW and WZ decay channels in which one W boson decays into an electron and a 
neutrino and the second W boson or the Z boson decays hadronically were studied in order to obtain an upper limit 
on the cross section and to restrict possible anomalies in the coupling parameters. In these processes, the W and Z 
bosons that decay hadronically to produce two jets in the detector cannot be differentiated due to the limitations of jet 
energy resolution. Figure || shows the Pt spectrum for W bosons in WW and WZ production from the leading-order 
theoretical prediction [ p5[ . This paper describes in detail an analysis summarized in Ref . [Q , in which the Et spectrum 
of the W bosons, produced with two or more jets which could have come from a hadronic W or Z boson decay, was 
compared to the expected SM signal plus background to set limits on anomalous WWj and WWZ couplings. The 
CDF Collaboration has studied the jet-jet decay mode and reported ||l^ limits on anomalous WW^ and WWZ 
couplings. 

A new result on the anomalous couplings from a combined fit is presented. Since the Wj, WW to dileptons, and 
WW/WZ to electron plus jets analyses measured the same coupling parameters, a combined fit to all three data sets 
was performed, yielding improved limits on the anomalous couplings compared to the individual analyses. 



B. ^77 and ZZy Trilinear Gauge Boson Couplings 



The interactions of pairs of neutral gauge bosons, the Z boson and the photon, can be studied through the Zj 
production process. The Feynman diagrams for the qq' — > Zj processes are obtained by substituting Vq = 7 or Z, 
Vi — Z and V2 = 7 in Figs. 0(a)-(c). There are no ZZj and .^77 couplings of the type shown in Fig. ^(c) in the 
SM; thus, there is no destructive interference of the u and t-channel amplitudes, such as occurs in Wj, WW, and 
WZ production. A Z boson is identified by its leptonic decay products, a pair of high px leptons {£ = e or fj,), or 
by the imbalance of momentum in the event due to not detecting the neutrino pair. The Drell-Yan production of 
a Z boson or virtual photon, followed by radiation of a photon off the charged Icpton from the Z boson or virtual 
photon decay products, also contributes to the charged lepton final states, as shown in Fig. ^. As with the final state 
radiation from W boson decay products, the photon from the Z boson decay products is preferentially emitted along 
the charged lepton direction; the process can be suppressed by imposing a cut on the separation between the charged 
lepton and the photon. The most general Lorentz and gauge invariant ZVj vertex is described by eight coupling 
parameters hY{i = 1, ...,4) ||l^. The anomalous couplings are parameterized as form factors hY — /iio/(l + s/A^)", 
where s is the square of the invariant mass of the Zj system, n = 3 for hY^, and n = 4 for hY 4- This is discussed in 
more detail in Appendix 1. Figure u] shows the Et spectra of photons predicted for the SM and the non-SM model 
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production processes. The distributions for the non-SM ZZ^ and Z77 couphngs exhibit a large increase of the cross 
section at high Et- Limits on the anomalous coupling parameters are obtained from a maximum likelihood fit to 
the Et spectrum of the photons, as in production. Previously, CDF has published limits on the ZZ^ and Z77 
anomalous couplings [is| ] using the ee and /i/i final states. Recent results from the D0 experiment are presented here 
in more detail than in the previous pjq| summaries and include the ee, /i/i, and neutrino final states. 



C. Outline of Paper 

The paper is organized as follows. In Section Othe D0 detector and the techniques used to identify particles from 



the collisions are discussed. Section [II is a summary of various Monte Carlo modelling tools used in these analyses. 
Section |^ discusses the 1992-1993 collider run and data samples. Section ^ describes a measurement of the WW^ 
coupling parameters using 14^7 events where the W boson decays into a high px electron or muon and a neutrino. In 
Section^ the results of a search of WW {iv){i'v') process are presented. This is followed by a description and the 



results of the analysis for WW and WZ production with subsequent decay to ev and at least two jets in Section VIl 
Section VIII| describes the combined limits on the anomalous couplings from all of the W^ and WW/WZ analyses 



Section IX presents a measurement of ZZ'y and ^77 coupling parameters using Zj production events where the 
Z boson decays into ee, fifi or vD. The results are summarized and reviewed in Section Appendix 1 contains 
a detailed review of trilinear gauge boson couplings and the effective Lagrangian. Finally, Appendix 2 contains a 
discussion of the binned likelihood fitting procedure used in determining the anomalous couplings limits. 



II. DETECTOR 

The major components of the D0 detector were a nonmagnetic central tracking detector system, a hermetic 
liquid-argon uranium calorimeter and a muon spectrometer system with a toroidal magnetic field. A perspective view 
of the detector is shown in Fig. depicting the three major systems. 

The central tracking detectors (CD), shown in Fig. ^, included the Vertex Drift Chamber (VTX), the Transition 
Radiation Detector (TRD), the Central Drift Chamber (CDC) and two Forward Drift Chambers (FDC's). The VTX, 
TRD and CDC were arranged in three concentric layers around the beamline, from the beampipe out to the central 
calorimeter. The wires in the FDC's were oriented perpendicular to the beamline. The entire CD was contained in the 
cylindrical volume (r — 76 cm, z — ±135 cm) bounded by the calorimeter cryostats. The CD measured the trajectory 
of charged particles with a resolution of 2.5 mrad in (f> and 28 mrad in 6, where (j) and 9 are the azimuthal and polar 
angles of the track, respectively, and covered the region |?7| < 3.2 in pseudorapidity, where 77 = — ln(tan|). From these 
measurements, the position of the interaction vertex along the beam direction was determined with a resolution of 8 
mm. The presence of a CD track or hits pointing towards a shower was the key element for distinguishing electrons 
from photons. The CDC and FDC's also provided ionization energy loss measurement for separating single electrons 
from closely-spaced photon conversion pairs where the photon converted before it reached the tracking detector. 

The calorimeter was divided into three parts, each enclosed in a steel cryostat to contain the liquid argon: the central 
calorimeter (CC) and two end calorimeters (EC's) as shown in Fig. Each consisted of an inner electromagnetic (EM) 
section, a finely-segmented hadronic section (FH) and a coarsely-segmented hadronic section (CH). The scintillator- 
based intercryostat detectors (ICD's), which improved the energy resolution for jets that straddled the central and 
end calorimeters, were inserted into the space between the cryostats. The calorimeter covered the pseudorapidity 
range \r]\ < 4.2. The boundaries between the CC and EC's were chosen to be approximately perpendicular to the 
beam direction and to match the transition between the CDC and FDC's. 

Each EM section was divided into four longitudinal layers forming a thickness of 21 radiation lengths. The hadronic 
sections were divided into four (CC) or five (EC) layers and were 7-10 nuclear interaction lengths thick. The calorime- 
ter was transversely segmented into projective towers with A77 x A(p = 0.1 x 0.1, where is the azimuthal angle. The 
third layer of the EM calorimeters, where the maximum energy deposition from EM showers was expected to occur, 
was segmented more finely into cells with A77 x A(f) ~ 0.05 x 0.05. The azimuthal position resolution for electrons 
with energy above 50 GeV was approximately 2.5 mm. 

The calorimeter provided the energy measurement for electrons, photons, charged hadrons, and jets. The energy 
resolution of the D0 calorimeter was measured in a test beam for electrons and pions. The energy resolution for 
electrons and photons was (Je/E — ^jE{GeV) © 0.4%. For charged pions and jets the resolutions were approx- 



imately 50%/v/£^(GeV) and E{GeY), respectively pIpTl. The transverse energy of a neutrino was inferred 



from the undetected transverse energy, which is the negative of the vector sum of all the transverse energies in 
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the calorimeter cells. Using a minimum bias data sample, the resolution for each component of the missing transverse 
energy, If^ and Ify, was measured to be 1.1 GeV + 0.02( ^ -Et), where ^ Et is the scalar sum of transverse energies in 
all calorimeter cells. For one analysis (see Section IX C| ) the I^t was calculated from the negative of the vector sum of 
the transverse energies in towers with Et > 200(400) MeV in the EM (FHl) calorimeters. The hadronic calorimeter 
scale was determined by comparing the transverse energy of the recoil against that of the electron pair in Z ee + X 

events. The resolution of the x and y components of the I^t was a = \J (4.88 GeV)^ + (1.34 • 

The muon spectrometer system, shown in Fig. H, consisted of solid-iron toroidal magnets and sets of proportional 
drift tubes (PDT's). It provided identification of muons and determination of their trajectories and momenta. Since 
muons from decays of W and Z bosons are primarily in the central region, the analyses presented here used only the 
wide angle muon spectrometer (WAMUS) which subtended the region \rj\ < 2.5. The WAMUS system consisted of 
three layers: the "A-layer" with four planes of PDT's, located between the calorimeter and the toroid magnets; and the 
"B-" and "C-layers" each with three planes of PDT's, located outside the toroid magnets. The toroids were magnetized 
to ±1.9 T. The wires in the drift tubes and the magnetic field lines in the toroids were oriented transversely to the 
beam direction. The muon system mounted on the central (forward) muon toroid, covering approximately |r/| < 1 
(1 < \r]\ < 1.7) is referred to as the "CF (EF)" region. 

The total material in the calorimeter and iron toroids varied from 13 to 19 interaction lengths, making background 
from hadronic punchthrough negligible. Because of the small tracking volume, the background to prompt muons from 
in-flight decays of tt and K mesons was also negligible. 

The muon momentum p was determined from its deflection angle in the magnetic field of the toroid. The momentum 
resolution was limited by multiple scattering in the calorimeter and toroid, knowledge of the magnetic field integral, 
and the accuracy of the deflection angle measurement. The momentum resolution, determined from J/iJj ^/z and 
Z fifj, events, was cr(l/p) = 0.18(p — 2)/p^ ® 0.008 {p in GeV/c), where © indicates addition in quadrature. 



A. Trigger 

A multi-level, multi-detector trigger system was used for selecting interesting events and recording them to tape. 
A coincidence between hits in two hodoscopes of scintillation counters (level 0), centered around the beampipe, was 
required in order to register the presence of an inelastic collision. These counters also served as the luminosity monitor 
for the experiment. The level 1 and level 1.5 triggers were programmable hardware triggers which made decisions 
based on combinations of detector-speciflc algorithms. The level 2 trigger was a farm of 48 Vax 4000/60 and 4000/90 
computers which filtered the events based on reconstruction of the information available from the front-end electronics. 
Losses from the Main Ring beam, usually involved in the production of antiprotons, caused backgrounds in the muon 
system and calorimeter. In the analyses presented here, triggers which occurred at the times when a Main Ring proton 
bunch passed through the detector were not used. Similarly, triggers which occurred during the first 0.4 seconds of 
the 2.4 second antiproton production cycle were vetoed. These "Main Ring vetoes" accounted for approximately 25% 
trigger deadtime. 

At each level of the trigger, the D0 trigger system gathered the results from each of the detector-specific triggers 
and filters. In this way, trigger decisions could be made from combinations of different detector-specific results. Table || 
is a compilation of the triggers used in the various analyses presented in this paper. 



1. Calorimeter Trigger 

The level 1 triggers for electromagnetic showers were based on analog sums of transverse energy in calorimeter towers 
with A77 X A(j) — 0.2 X 0.2 and with two longitudinal sections, EM and FH. The level 1 EM trigger required transverse 
energy in the EM section of a trigger tower to be above preselected thresholds. The level 2 EM algorithm identified 
electrons and photons by forming clusters, around level 1 trigger towers, of transverse energy read out from the four 
layers of the EM calorimeter and the first layer of the FH calorimeter. The clusters were of size A77 x A(p — 0.3 x 0.3, 
centered on the highest Et tower in the cluster. The longitudinal and transverse profile of the cluster had to satisfy 
the following requirements which were designed to discriminate electrons and photons from hadronic showers. The 
fraction of the cluster energy in the EM section had to exceed a value which depended on the energy and location of 
the cluster in the calorimeter. The transverse shape classification was based on the energy deposition pattern in the 
third EM layer. The difference of the energy depositions in two regions, of size Ar/ x = 0.25 x 0.25 and 0.15 x 0.15 
and centered around the cell with highest Et, had to be within a range which depended on the total cluster energy. 
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Another calorimeter based trigger was the "Missing Et" trigger. At level 2 the I^t was formed from the negative 
of the vector sum of the Et deposited in the calorimeter and ICD cells, corrected for the vertex position. 

2. Muon Trigger 

The muon level 1 and level 1.5 triggers required coincidences of hits in the PDT's consistent with a muon originating 
from the collision region. The level 1 algorithm combined coincidences of hits in PDT cells into 60 cm- wide hodoscopic 
elements. If a combination of hodoscopic elements matched a preprogrammed pattern of a muon track, the event 
was accepted. In the central region, three layers of PDT's, each with at least two hit planes, were required except in 
regions where detector services and support limited the coverage of the one of the layers. In the forward region, defined 
approximately as 1.0 < |?7| < 2.5, three layers of PDT's were required, with at least three hit planes in the A-layer and 
two hit planes in both the B- and C-layers. The trigger required a minimum pT of 3 GeV/c and became fully efficient 
at about 6 GeV/c. The level 1 trigger efficiency was (79 ± 3)% for the central region and (36 ± 12)% for the forward 
region. At level 1.5, the hodoscopic elements had half-cell resolution, providing a sharper px turn-on. A three-layer 
requirement made at this trigger level reduced the acceptance of the central muon system by approximately 15%. 

At level 2, the first stage of the muon reconstruction algorithm, which consisted of the pattern recognition and 
initial track fit, was performed. To minimize processing time, the search for muon candidates in the forward region 
was restricted to the sectors which had a level 1 trigger. A valid level 2 trigger was a three-dimensional muon 
track with hits in at least two planes of two PDT's. The level 2 muon trigger program calculated several quantities 
that provided information on the quality of the muon track including: the goodness of track fit in the PDT drift 
view and along the PDT wire, the projections of the track to the interaction point in both views and the number 
of hits used to fit the track. A track quality variable was defined as the number of these quantities that failed 
the standard criteria. In addition, in the forward region, muon candidates formed with less than six hits on the 
track were discarded, since they were likely to be random background hits in coincidence. In the central region, 
cosmic ray muons were identified if there was evidence of a single muon penetrating the entire detector; muon 
candidates with a track within 20 deg in and 5 deg in or hits within 60 cm (roughly 5 deg) of the projection 
of the muon track into the opposite side PDT's were rejected. The muon was accepted by level 2 if the pt was 
above the desired threshold and if the track quality variable was zero ("tight" standards), or one ("loose" muon 
standards). The muon level 2 trigger efficiency was determined to be (95 ± 3)% excluding effects of the chamber 
efficiencies and geometrical acceptance for the "loose" muon requirements of the MU-ELE trigger (see Table |). 



Trigger Name 


Level 1 


Level 1.5 


Level 2 


Analyses 


MU-MAX 


1 /i, \ri\ < 1.7 


1 \v\ < 1.7 


1 (tight), PT > 15 GeV 


WW fi^i 


MU-ELE 


1 EM tower, Et > 7 GeV 
1 M, < 1.7 




1 e or 7, £t > 7 GeV 
1 jj. (loose), Pt > 5 GeV/c 


WW efi 
Wj — > nv-y 
Z'y fi^i'Y 


ELE-HIGH 


1 EM tower, Et > 14 GeV 




1 e or 7, £t > 20 GeV 


WW/WZ 
ev jet jet 
Z7 — > vv^ 


ELE-MAX 


1 EM tower, Et > 10 GeV 




1 e or 7, £t > 20 GeV 
^T > 20 GeV 


W"f ev^ 


ELE-2-HIGH 


2 EM towers, Et > 7 GeV 




2 c and/or 7, Et > 10 GeV 


WW ee 


ELE-2-MAX 


2 EM towers, Et > 7 GeV 




2 e and/or 7, Et > 20 GeV 


Z') — > 667 



TABLE I. Triggers used in the analyses presented in this paper. The ELE-2-MAX trigger was a subset of the ELE-2-HIGH 
trigger which included shower shape cuts on the EM candidates. 



B. Muon Identification 

Muons were identified as tracks in the muon PDT's in association with tracks in the CD and energy deposits in the 
calorimeter. The momentum of the muon was computed from the deflection of the track in the magnetized toroid. The 
track fit used a least-squares calculation which considered seven parameters: four describing the position and angle 
of the track before the calorimeter, two describing the effects due to multiple Coulomb scattering, and the inverse of 
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the muon momentum, \/p. This seven parameter fit was applied to 16 data points: vertex position measurements in 
the X and y directions, the angles and positions of the track segments before and after the calorimeter and outside 
of the toroid magnet, and two angles representing the multiple scattering of the muon in the calorimeter. The fit 
determined the charge of the muon and which CD track, if any, matched the muon. The muon momentum was then 
corrected for the energy lost in the calorimeter using a GEANT-based |2^] detector model. 

In the following, the quantities used to describe the muon tracks are presented. The definitions for muons differ 
slightly among the various analyses because of the nature and magnitude of the backgrounds. Table || lists the five 
different definitions of muons in the analyses described in this paper. 

1. Muon Track Quality 

The muon reconstruction algorithm defined a muon track quality, similar to that used in the level 2 trigger, which 
contained information about the number of hits on the track from each layer of muon PDT's, the track impact 
parameters, and the goodness of the track fit. If the track did not satisfy criteria on more than one of the above 
quantities, the muon candidate was rejected. Figures |ll|(a) and (b) show the impact parameters in the track bend 
view (r — z plane), tbcnd, and in the track nonbend view {x — y plane), fononbcnd, for muons which satisfied all of the 
other selection criteria. The 3D impact parameter was sometimes used in lieu of the combination of the r — z and 
X — y selection criteria. 

2. Fiducial Requirements 

Muons which passed through the region between the CF and EF toroid magnets near \'q\ « 0.9 may have traversed 
a smaller amount of magnetized iron and thus have a reduced momentum resolution. To reject these poorly measured 
muons, all of the muon identification definitions in this paper except "Loose 11" required the minimum magnetic field 
integral along the muon track, J Bdl, to be at least 2.0 Tesla-meters. 

The "Tight 11" , "Tight III" , and "Loose 11" definitions required that the muon have hits in the A-laycr, between 
the calorimeter and the toroid magnet. Making this requirement reduced the fake- muon background in the forward 
region but also reduced the acceptance by limiting the pseudorapidity coverage to approximately \ri\ < 1.7. A cut on 
\ri\ < 1.7 was used to restrict the muon tracks to those totally contained within the WAMUS spectrometer. 

3. Central Detector Track-Match 

Muon candidates were required to have a confirming track in the CD within a range in both the polar and azimuthal 
angles. This reduced the backgrounds from cosmic ray muons and from combinations of random hits. 

4- Calorimeter Confirmation 

Muons deposited energy in the calorimeter as they passed through it. It was required that at least 1 GeV of energy 
was deposited in the tower which contained the projection of the muon track through the calorimeter and the nearest 
neighboring towers. Figure pT|(c) shows the energy deposition in these towers around muons which passed the other 
selection requirements. 

5. Cosmic Ray Identification 

The muon track was refitted with the timing of the muon track with respect to the pp interaction as a floating 
parameter, tQ. This allowed identification of cosmic ray muons. Figure [TT|(d) shows the value of tg which resulted in 
the best track-fit for muons which passed all other selection requirements. The "Tight I" and "Tight 11" definitions 
required < 100 ns. 
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6. Muon Isolation 



Muons from the decay of pions, kaons, and heavy quarks were reduced by requiring that the muon be isolated 
from other jet activity. This was done in three ways. One isolation variable (2NN) was defined by summing the 
energy deposited in the calorimeter cells hit by the muon and two nearest neighbors, subtracting the energy expected 
to have been deposited by the muon, and dividing the difference by the uncertainty. This was required to be less 
than five standard deviations. Another isolation variable (Halo) was defined as the difference between the energy 
deposited in a cone of size AR = 0.6 and the energy deposited in a cone of size AR = 0.2, where AR = ^J Arf + A(j)'^, 
around the muon in the calorimeter. This was required to be less than 8 GeV. The third isolation criterion (A_R^_jct 
) was that muons were spatially separated from the axis of any jet with Et > 10 GeV by at least AR — 0.5. 



Selection 


Tight I 


Tight II 


Tight III 


Loose I 


Loose II 


Analysis 


WW fifiuu 






WW efiuu 


Zj -> fifi-y 


Muon Quality 


V 


V 


V 


V 


V 


Back-to-Back 












Muons Removed 


V 




V 


V 


V 


Minimum 












Field Integral (T-m) 


2.0 


2.0 


1.9 


2.0 




3-Layers Required 






V 






A-Layer Required 




V 


V 




V 


Isolation 


2NN Cut 




A-R/j-jct 


A-Rfj-jct 


A-Rfi-jot 


Requirement 


Halo Cut 


> 0.5 


> 0.5 


> 0.5 


> 0.5 


Impact Parameter 


|3D1 < 22 cm 


\RZ\ < 22 cm 
\XY\ < 15 cm 


\RZ\ < 22 cm 


13D1 < 22cm 


\RZ\< 22cm 


to 


< 100 ns 


< 100 ns 








CD Match 


S(l} < 0.45 


5(f> < 0.25 


S(j} < 0.25 


S(j} < 0.45 


5(j> < 0.25 




561 < 0.45 


se < 0.30 


se < 0.30 


se < 0.45 


59 < 0.30 


Cal Confirm 


V 


V 


V 


V 


V 



TABLE II. Summary of the various muon identification definitions used in the analyses presented in this paper. 



C. Electron and Photon Identification 



Electrons and photons were identified by the properties of the shower in the calorimeter and the presence, or lack 
thereof, of a matching track in the CD. Using a nearest-neighbor algorithm, clusters were formed from adjacent EM 
towers containing significant energy deposition. The clusters for which the energy in the EM and first FH section 
of the calorimeter divided by the sum of the energies in the EM and all hadronic sections (EMF) was greater than 
or equal to 0.9 were flagged as possible electrons or photons. Figure |l^(a) shows the fraction of electrons from Z 
boson decays for which the EMF is above the value given on the abscissa. More detailed analysis of the calorimeter 
and tracking chamber information was then used to refine the sample as is described below. A summary of electron 



definitions is presented in Table III. The photon definitions are summarized in Table EV 



1 . Fiducial Coverage of EM Calorimeter 



All of the analyses presented in this paper made identical selection on the fiducial coverage of the EM calorimeter. It 
was required that the electron or photon have pseudorapidity within the range ±2.5. Furthermore, the EM calorimeter 
had a gap in the coverage at the transition between the CC and the EC. The four longitudinal layer coverage of the 
CC ended at pseudorapidity of 1.1 and the four longitudinal layer coverage of the EC returned at pseudorapidity 
of 1.5. Therefore, all analyses of photons and electrons in this paper made a fiducial selection which removed this 
transition region. The "Tight" photon identification criteria went one step further, requiring that photons in the CC 
have pseudorapidity within the range ±1.0. The samples used in the discussion of photon and electron identification 
presented below have these fiducial selections already applied. 
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2. Covariance Matrix 



The electron or photon shower shape was characterized by the fraction ol the cluster energy deposited into each 
layer of the calorimeter. These fractions were correlated, depending on the depth of the start of the shower and on the 
energy of the incident particle. In order to reject background using the shower shape, including these correlations, a 
comparison was made between the candidate and a reference sample of Monte Carlo electrons with energies ranging 
from 10 to 150 GeV. This comparison (H-matrix x^) was carried out in 41 observables: the fractional energies in 
layers 1, 2, and 4 of the EM calorimeter; the fractional energy in each cell of a 6 x 6 array of cells in layer 3 centered 
on the most energetic tower in the EM cluster; the logarithm of the total energy of the electron cluster, taking into 
account the depth dependence on the incident energy; and the position of the event vertex along the beam direction, 
taking into account the dependence of the shower shape on the incident angle. A separate reference shower shape was 
available as a function of 77, assuming (j) symmetry. Figure p^b) shows the fraction of electrons from Z boson decays 
for which the value of the H-matrix is less than the value given on the abscissa. Requiring that the H-matrix 
< 100(200) in the CC (EC) gave an efficiency of 94.9 ± 0.8% (100.01?;^%) for electrons with Et > 25 GeV. 

The efficiency for the H-matrix selection decreased if the Et fell below 25 GeV. The efficiency as a function of photon 
Et was measured in a test beam for both the CC and EC. This dependence was a dominant source of systematic 
uncertainty in the efficiency for low pt photons. Figure |l^ shows the efficiency versus pt for the H-matrix selection 
criteria for low pT "Loose" photons. 



3. Cluster Isolation 



The EM clusters were required to be isolated from other particles in the event in order to reduce the background 
from hadronic jets with high EM content. The isolation variable was 

_ i;(0.4) ~ EM{0.2) 

EM{0.2) ' ^ > 

where i?(0.4) was the energy deposited in all the calorimeter cells in a cone of radius R — 0.4 around the electron or 
photon and EM{0.2) was the energy deposited in the EM calorimeter in a cone of radius R — 0.2. For EM objects 
with Et < 20 GeV, there was deterioration of the efficiency of the fiso selection criteria. This was modeled with a 
turn-on curve in a way similar to the H-matrix efficiency described above. Figure p^(c) shows the fraction of electrons 
from Z boson decays passing an fiso selection criterion. Requiring fjso < 0.10 was 97.6 ± 0.6% (98.5 ± 1.4%) efficient 
for CC (EC) electron candidates. 



4. Electromagnetic Fraction 

The "Tight" photon identification criteria included the requirement that the energy deposited in the four EM layers 
be at least 96% of total energy in the calorimeter in a cone around the shower maximum. This was in addition to the 
EMF requirement discussed above. 



5. Electron Track Match 



Electrons and photons were distinguished from each other by the presence of a track consistent with the passage 
of a charged particle in the CD which pointed to the EM cluster in the calorimeter. An electron had such a track; 
a photon did not. The efficiency for track-finding was 86.7 ± 1.4% in the CDC and 86.1 ± 1.8% in the FDC's. By 
demanding a good spatial match between the cluster and the track, backgrounds due to accidental overlaps of charged 
particles with photons in EM jets were reduced. In the calorimeter, the shower centroid, Xc, was determined from the 
weighted mean of the coordinates xi of all cells containing the shower, Xc = ''^i^i/ ''^i- The weights were defined 
as Wi — max(0, + In^), where Ei was the energy in the ith cell, E the energy of the cluster, and wq a parameter 
chosen to optimize the position resolution. The logarithmic weighting was motivated by the exponential lateral profile 
of an electromagnetic shower. The azimuthal position matching resolution in the CC and EC was measured to be 
« 2.5 mm. The CD track was extrapolated to the shower centroid and the significance of the track match, TAI S, was 
formed between the position of the track and the centroid. For the CC this quantity was 
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TMScc = ^|i^r + i^r, (2.2) 



where was the azimuthal mismatch, Az the mismatch along the beam direction, and Sx was the resolution for 



the observable x. For the EC, Az was replaced by Ar, the mismatch transverse to the beam. Figure 12(d) shows the 
fraction of electrons from Z boson decays for which the track match significance variable is less than the criterion on 
the abscissa. Requiring TMS < 10 was 98.0 ± 0.6% (91.5 ± 1.8%) efiicient for CC (EC) electron candidates. 



6. Electron Track Ionization 



The tracks from e+e~ pairs produced in photon conversions due to interactions with material in the tracking 
chambers were often reconstructed as a single track. For such pairs, the ionization in the tracking chambers was 
expected to be twice that of a single charged particle. The distribution of ionization per unit length (dE/dx) for 
electrons from Z ^ ee decays and from EM clusters in an inclusive jet sample are shown in Figs. |l^(a) and p^(b). 
Most electrons had dE/dx « 1. The ionization in the inclusive jet sample shows a two-peaked structure. The lower 
peak, at dE/dx « 1 was due to single charged particles. The higher peak came from unresolved e+e" pairs. This 
background was rejected by removing electron candidates with dE/dx ^ 2. The veto requirement for CC(EC) was 
1.6 < dE/dx < 3.0 (1.6 < dE/dx < 2.6) and was 94.4 ± 1.1 (75.2 ± 3.7)% efficient for electrons from Z bosons. 



7. Loss of Photons due to Track Overlaps 



Some photons were mislabeled as electrons because of spatial overlap of the photon with a random track. The 
inefficiency introduced was estimated by looking for a track or tracks in a cone randomly oriented in but at the 
same rj location as the electrons in Z boson decays to ee. The assumption is that the probability of finding such a 
track or tracks at a given 77 is the same in Z boson production and double vector boson production. The probability 
for a random track overlap was found to be 6 ± 1% and 15 ± 1% for the CC and EC, respectively, the latter being 
higher due to the higher density of tracks in the forward direction. 



8. Photon Conversions 



Some photons were lost when they converted to e+e~ pairs in material in front of the CDC or FDC. The conversion 
probability was calculated using the GEANT simulation of the D0 detector. This probability depended on the 
pseudorapidity of the photon and is shown in Fig. |l5|. Averaged over the CC (EC) it amounted to a 10% (26%) loss 
of photons. There is a systematic uncertainty of 5%. 



9. Photon-Vertex Projection 



An algorithm, EMVTX, was developed to reduce the background from cosmic ray or beam-related muon 
bremsstrahlung which produced photons inconsistent with having originated at the event vertex. The energy- weighted 
centers of the cluster in each of the four layers of the EM calorimeter plus the vertex position, and their uncertainties, 
were used to compute two dimensional fits. The resulting was then converted into a probability for the photon to 
have originated at the vertex. It was required that the probability of the RZ and XY projections, Puz and Pxy, each 
exceed 1%. Comparison of the Prz and Pxy distributions for electrons from Z bosons and from photons resulting 
from cosmic ray bremsstrahlung are shown in Fig. |l6[ In case there were multiple vertices in the event, the one with 
the highest Puz was selected as the vertex for the interaction. This vertex was then used in computing the missing 
transverse energy and E'^. The vertex resolution provided by this algorithm was approximately 17 (11) cm in the RZ 
{XY) planes. 
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10. Hits Along Photon Roads 



The backgrounds to photons from electrons and high-EM content jets with unreconstructed tracks were reduced by 
looking for hits in narrow roads between the vertex and the EM cluster in the calorimeter. In particular, a background 
to Zj i/i/j was W ^ ev events where the electron was misidentified as a photon. The tracking algorithm could 
have been confused by extra hits or have missed the track because the wrong vertex was selected (in the case of 
multiple vertices) or because the vertex was reconstructed poorly and the track pointed in a different direction. An 
algorithm was developed, called HITSINFO, to identify this background. Roads were defined between the cluster and 
each reconstructed vertex as well as the vertex position obtained by the EMVTX algorithm. The road size depended 
on the tracking chamber. The following road sizes were used: 

MvTX = 0.005, AcpvTX = 0.012; 
A6cDC = 0.050, A^cDC = 0.0075; 
A0FDC = 0.005, A(j)FDC = 0.015; 

where the angles, and Acp, were the half-opening angles of the roads in the RZ and XY planes. 

The roads were examined for tracks and hits. The photon candidate was required to have no tracks from 
any vertex. Further requirements were made on the fraction of available wires hit, the number of reconstructed 
track segments, and on the number of hits, depending on the tracking sub-detector. The selection criteria were 
optimized using the Z ^ ee sample with one of the electrons being misidentified as a photon due to track- 
ing chamber ineflaciency. The efficiency was calculated using a sample of "emulated" photons obtained by rota- 
tion of the positions of the electron energy clusters by 90 degrees in cp and then applying the selection criteria. 



Selection 


Tight I 


Tight II 


Tight III 


Tight IV 


Loose I 


Analysis 


WW ee and e/i 


Z7 667 


VK7 — > ewy 


WW/WZ ev jet jet 


Z7 667 


EM Fraction 


> 0.90 


> 0.90 


> 0.90 


> 0.90 


> 0.90 


Track Match 












CC (EC) 


< 100(100) 


< 100(200) 


< 100(200) 


< 100(100) 


< 100(200) 


Isolation 


< 0.1 


< 0.1 


< 0.15 


< 0.1 


< 0.1 


TMS 


< 10(7 


< 10(7 


< 10(7 


< 5(7 




dE/di- 


/ 

V 










CC Efficiency 


72.9 ± 2.3% 


78.1 ±2.3% 


79 ± 2% 


76.7 ± 1.6% 


90.2 ± 1.3% 


EC Efficiency 


51.0 ±3.6% 


70.8 ± 3.4% 


78 ± 3% 


62.0 ±3.1% 


97.1 ±2.9% 



TABLE III. Summary of the various electron identification definitions used in the analyses presented in this paper. 



Selection 


Tight 


Loose 


Analysis 


Z^t lyiy^i 


Wj 

Z'l (..(.-, and ///i7 


EM Fraction 


> 0.96 


> 0.90 


CC (EC) 


< 100(100) 


< 100(200) 


Isolation 


< 0.1 


< 0.1 


Matching Track 


Veto 


Veto 


EMVTX 


V 




HITSINFO 


V 




CC Efficiency 


57 ± 2% 


74 ± 7% 


EC Efficiency 


56 ± 4% 


58 ± 5% 



TABLE IV. Summary of the various photon identification definitions used in the analyses presented in this paper. The 
"Loose" efficiencies do not include the pt dependent effects important at low Pt- Similarly, the "Tight" efficiency applies to 
only the high pT photons within the fiducial region. 
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D. Jet Reconstruction 



Jets were reconstructed using cone algorithms with cone sizes, ATZ, of 0.3, 0.5, and 0.7 for the analyses presented 
in this paper. The algorithm was as follows. First, a jet candidate was identified by forming prcclustcrs of size 
At] X Ac6 = 0.3 X 0.3, centered on the highest Et tower in the cluster, from a list of jet towers with Et > 1.0 GeV 
ordered by Et- Next, the jet direction was determined by an iterative process. A cone of size ATZ was placed around 
a new Et wcightc^d jet c;c!nt(;r of towers and the process was repeated until the jet direction became stable. If two jets 
shared energy, they were combined or split, based on the fraction of energy shared relative to the Et of lower Et jet. 
If the shared energy was greater than 50% of the lower Et jet, the jets were merged. 

The; jet ciiiergy was corrected for a number of effects. These included energy contributed to the jet from the 
underlying event, energy from the jet which escaped the jet cone, energy lost due to the zero-suppression, as well as 
the overall jet energy scale. A cone of radius ATZ = 0.7 was selected by all analyses presented here, except for the 
WW/WZ ev jet jet analysis, where a cone of radius ATZ = 0.3 was used, and the Zj vv^f analysis, where a cone 
of radius ATZ = 0.5 was used. The small cone size was advantageous for detecting two closely-spaced jets expected 
from high-pr W boson decays. The larger cone size had smaller Et corrections. 

E. Missing Transverse Energy 

The missing transverse energy in the calorimeter If,^"^ was defined as 

^cal = ^^cal^+^cal^, (2.3) 

where 

= -Y.E, sin cos 0i - ^ AEl, (2.4) 

i 3 

9tI = - E ^^"^ ' (2-5) 

i 3 

The first sum is over all the cells in the calorimeter and ICD. The second sum is over all the corrections in Et applied 
to all the electrons and jets in the event. In order to obtain the best resolution, the corrections AElp were those from 
reconstructing the event with a jet of cone size ATZ = 0.7. 

The soiirces of I^t included neutrinos, which escaped undetected, and the energy imbalance due to the resolution 
of the calorimeter and muon system. The missing transverse energy was corrected if there were muons in the event. 
The transverse momenta of the muons were removed from the to form the total missing Et, whose components 
were: 

VT. = rTi-J2P^.% (2.6) 



^t, = ^^^^-EK'- (2-7) 

i 

In that which follows in this paper, analyses not involving muons did not distinguish between ^j?' and I^t', the "cal" 
superscript is then ignored. 

III. MONTE CARLO SIMULATION 

In order to determine effects due to experimental limitations such as detector acceptances, resolutions, and efficien- 
cies on the expected signal and background, and to provide a cross-check for many of the quantities measured with 
the data, simulations of the detector and trigger were developed. Various levels of sophistication were used, depending 
on the detail required. 
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A. Detector simulation programs 



The most detailed model of the detector was the GEANT |22| simulation. The D0 implementation of GEANT, 
D0GEANT, included details of the geometry of individual detectors, instrumental efficiencies and resolutions, and 
particle responses. The performance of the D0GEANT program was confirmed by comparing the simulation results 
with the data taken from test beams [pT| , cosmic ray muons and pp collisions. It was typically used to predict and 
cross-check the effect of variations in the particle identification requirements on the efficiency for leptons and jets. 
It was also used to predict and cross-check the effect of changing the kinematic requirements on the number and 
characteristics of some of the signals and backgrounds. 

In a typical application, an event generator such as PYTHIA |23) or ISAJET was used to create a list of particles 
produced in the collision. The simulation converted this into a Monte Carlo event with the same format as the 
digitized information from the real collision. This Monte Carlo event was then reconstructed in the same way as the 
data. As it was a very detailed detector simulation, it consumed relatively large quantities of computer resources. 
This limited its application to problems of manageable scale. 

In order to speed-up the GEANT simulation, the calorimeter response for electrons, photons, and hadrons could 
be modeled using a database of particle showers called the Shower Library. The Shower Library ^ was created by 
storing the energy deposition in each calorimeter cell for each shower that was generated using GEANT in the full 
shower mode. Each shower was stored in a list together with its particle identity, momentum, pseudorapidity, azimuth, 
and collision vertex origin. When using the Shower Library to simulate the response of a particle in a Monte Carlo 
event, a shower of the appropriate type was selected randomly from the library and added to the event. This method 
was use ful, for example, in determining the efficiency for dijets in the WW/WZ — > ev jet jet analysis presented in 
Section VII, where the advantage of speed made it possible to create a parameterization of the efficiency. 

An even faster simulation, D0FAST, with correspondingly less detail, used simplified geometrical structures of the 
D0 detector and parameterizations of the detector response including energy (momentum) resolutions, particle iden- 
tification efficiencies, and trigger turn-on curves obtained from the data and described in the previous section. Careful 
comparisons were performed between D0FAST and D0GEANT for the processes with Standard Model couplings to 
ensure that nothing important was lost in using the former. This simulation was used, for example, to model the 
acceptance for the grid points in the anomalous coupling parameter space, where dozens of grid points were used, 
each with 10,000 to 100,000 MC events. A similar fast Monte Carlo was used to estimate the background from Z 
boson decays for the WW dileptons analyses. 



B. Trigger Simulation 



In order to optimize and cross-check the efficiency of the triggers described in Table |, and to provide a method for 
finding the efficiency of combinations of separate triggers, a detailed simulation of the trigger algorithm was made. The 
list of available triggers, particularly those used for monitoring the higher-p^ triggers, changed from time-to-time over 
the course of the run as the luminosity increased. Occasionally the algorithms were improved as our understanding of 
the detector improved. The trigger simulation, TRIGSIM, was used to pre-test the changes in the trigger. The output 
from the GEANT simulation was processed by the simulator using level 1 (Ll.O) and level 1.5 (LI. 5) hardware and 
level 2 (L2.0) software simulations. The L2.0 simulation used software identical to that used in the L2 computers. 
The results were then compared to the arrays of available triggers and the events were marked as passed or failed. 
The simulator was cross-checked against the actual trigger using real data as input events. 



IV. DATA SAMPLES 



During the 1992-1993 pp collider run, the Fermilab Tevatron, operating at a center of mass energy of y/s = 1.8 
TeV, delivered a total integrated luminosity of / Ldt = 21.8 pb~^. Typical instantaneous luminosities of 4 x 10'^° 
cm~^sec~^ were attained. D0 collected 14.4 pb~^ to tape. The difference between delivered and collected luminosity 
was dominated by the dead time incurred due to operation of the Main Ring accelerator. A small part of the data 
was lost due to operational difficulties and hardware problems (bad runs) at the time of data collection. 

The luminosity was calculated by measuring the rate for pp nondiffractive inelastic collisions using the level 
scintillation counter hodoscopes. The normalization for the luminosity measurement and the 5.4% systematic un- 
certainty came from the pp inelastic cross section and the uncertainty in the acceptance of the counters p6|. The 
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final integrated luminosity varied from trigger to trigger for a number of reasons. The muon LI. 5 triggers started 
operating approximately six weeks after the muon Ll.O and calorimeter triggers. The muon triggers tended to 
be prescaled at high luminosities because they had higher Ll.O and LI. 5 trigger rates than the calorimeter trig- 
gers. Finally the analyses which used only EM objects could use luminosity collected while the muon system 
had hardware problems whereas the muon system, which relied on the calorimeter as part of muon identifica- 
tion, could not use luminosity collected when the calorimeter had a problem. The luminosity for a given trig- 
ger may have varied slightly from analysis to analysis depending on the bad run list used. Tabic ^ shows the 
total integrated luminosity, after bad run removal, for each trigger used in the analyses presented in this paper. 



Trigger 


/ Ldt (ph-^) 


MU-MAX 


12.2 ±0.7 


MU-ELE 


13.8 ±0.7 


ELE-HIGH 


13.7 ±0.7 (13.1 ±0.7) 


ELE-MAX 


13.8 ±0.7 


ELE-2-HIGH 


14.3 ±0.8 


ELE-2-MAX 


14.3 ±0.8 



TABLE V. Integrated luminosity for each trigger after accounting for the effects of the Main Ring and for bad runs due to 
hardware problems. The ELE-HIGH trigger has separate luminosities depending on whether the calorimeter or calorimeter 
plus muon system were checked. 



V. Wy ANALYSIS 



A measurement of the WWj couplings using pp —>■ tv^ + X {£ — e, fi) events is presented in this section. These 
events contained the production processes, pp W-f, followed hy W ^ or the final state radiation process 
W ^ ^ iv^, as shown in Figs. ^ and ^. Anomalous coupling parameters would enhance the production cross 
section, leading to an excess of events with high transverse energy photons, well separated from the charged lepton. 
Figure || shows the photon spectrum from SM and anomalous couplings predicted by theory |jl^ . The procedure of 
the analysis was to obtain a candidate sample, estimate the background contribution as a function of photon Et, and 
compare the background-subtracted candidate photon spectrum with that expected from various anomalous WWj 
couplings. In the following, the electron and muon channels are referred to as W{ev)^ and W{^v)^^ respectively. 



A. Event selection 



The W{ev)^ candidates were obtained by searching for events with an isolated high Et electron, large missing 
transverse energy, and an isolated photon. The data sample was taken with a Ll.O trigger that required at least 
one EM tower with Et > 10 GeV and the ELE-MAX trigger at level 2, that required an isolated EM cluster with 
Et > 20 GeV and I^t > 20 GeV, as described in Table |. The data sample corresponded to an integrated luminosity 



of 13.8 ± 0.7 pb . The electron was required to pass the "Tight III" requirements of Table III and the photon to 



pass the "Loose" requirements as described in Table IV. The electron and photon were required to be within the 



fiducial region of the calorimeter, as discussed in Section II CI, and at least 0.01 radians away from the azimuthal 



boundaries of the 32 EM modules in the CC. Kinematic selection was made requiring E^ > 25 GeV, I^t > 25 GeV, 
and Mt > 40 GeV/c^, where Mt is the transverse mass of the electron and I^t vector defined as 

Mr = [2E^$Til - cos^^")]^/^ , (5.1) 

and (p^" is the angle between the electron Et and the I^t- 

The W{fii')^ candidates were obtained by searching for events with an isolated high pT muon and an isolated 
photon in the data sample taken with the MU-ELE trigger described in Table ||. The sample corresponded to an 
integrated luminosity of 13.8 ± 0.7 pb^^. A muon track satisfying the "Tight 11" definition of Table |l| was required. 
Kinematic selection was made requiring > 15 GeV/c and I^t > 15 GeV. To reduce background from Zj events, 
where the I^t resulted from muon pT mismeasurement, events were rejected if they contained an additional muon 
track with > 8 GeV/c. 
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The requirements on photons were the same for both the electron and muon samples. The photon was required 
to have > 10 GeV. The separation between the photon and charged lepton, A7?.^-y, was required to be > 0.7. 
This requirement suppressed the contribution of the final state radiation process, and minimized the probability for a 
photon cluster to merge with a nearby calorimeter cluster associated with an electron or a muon. The above selection 
criteria yielded 11 W{ei')j candidates and 12 W(pv)j candidates. 



B. Efficiencies 



The trigger and offline lepton selection efficiencies, shown in Table [VI| , were primarily determined using Z ^ ££ 
events, requiring only one of the leptons to pass the trigger and selection criteria. Thus the second lepton provided an 
unbiased sample to measure efficiencies. The efficiency for the I^t requirement of the ELE-MAX trigger was calculated 
using the events which passed the ELE-HIGH trigger, which had no I^t requirement. The detection efficiency of the 
photons with Et > 25 GeV was determined using electrons from Z decays. For photons with lower Et there was a 
decrease in detection efficiency due to the cluster shape requirement, determined using test beam electrons, and the 
isolation requirement, which was determined by measuring the energy in a cone of radius TZ = 0.4 rotated randomly 
in azimuth in the inclusive W{eu) sample. Combining this i?T~dependent efficiency with the probability of losing 
a photon due to e+e~ pair conversion, 0.10 (0.26) in the CC (EC), and due to an overlap with a random track 
in the event, with probability 0.065 (0.155), the overall photon selection efficiency was estimated to be 0.43 ± 0.04 
(0.38 ± 0.03) at £:J = 10 GeV which increased to 0.74 ± 0.07 (0.58 ± 0.05) for > 25 GeV. 

The kinematic and geometrical acceptance was calculated as a function of coupling parameters using the Monte 



Carlo program of Baur and Zeppenfeld |10|, in which the 14^7 production and radiative decay processes were gen- 
erated to leading order, and higher order QCD effects were approximated by a K-factor of 1.34. The MRSD— ' 
structure functions were used and the px distribution of the system was simulated using the observed pr 



spectrum of the W in the inclusive W{ev) sample. Table Vl lists the acceptances for the SM production of W{iv)^ 





\v\ < 1-1 


M^(ei')7 
Erie) > 25 GeV 

1.5 < \v\ < 2.5 


\v\ < 1.0 


W(/ii^)7 
prifi) >15 GeV/c 

1.0 < \rj\ < 1.7 


e« 

SM 
^Acc 


0.98 ±0.02 
0.79 ± 0.02 


0.98 ±0.02 
0.78 ± 0.03 

0.11 ±0.01 


0.74 ±0.06 
0.54 ±0.04 


0.35 ±0.14 
0.22 ±0.07 

0.29 ± 0.02 



TABLE VI. Summary of trigger (etrig) and lepton selection (e^) efficiencies and geometrical acceptances (e^^c) foi" the SM 
Wy production events. 



C. Backgrounds 



The background estimate, summarized in Table VII, included contributions from: Z7, where the Z decays to 
££, and one of the leptons was undetected or was mismeasured by the detector and contributed to with 
W ^ TV followed by T ^ ivv] and W + jet(s), where a jet was misidentified as a photon. The backgrounds due to 
Z7 were estimated using the Zy event generator of Baur and Berger followed by a full detector simulation using 
the GEANT program j^. It should be noted that a{Z{U)y)la{W{(.v)-f) is about 0.5 (rather than 0.1 which is the 
ratio of cross sections oi Z ^ tl and W tv), since the W{£v)y process is suppressed by interference between the 
production diagrams and since the Z boson has twice as many leptons from which a photon can be radiated. The 
background due to Wy {tv)^ was estimated from the ratio of the detection efficiencies oi W ^ tv ^ evv and 
W ^ eu processes. The ratio was found to be 0.019 ± 0.002, using the ISAJET [|2^ event generator followed by the 
GEANT detector simulation. 

The W + jets background was estimated using the probability, 7'(j — > "7"), for a jet to be misidentified as a 
photon. The probability was determined as a function of Et of the jet by measuring the fraction of jets in a sample 
of multi-jet events that passed the photon identification requirements. Of course, some of the "fake rate" was due to 
real photons in the jet sample. The fraction of direct photon events in the multi-jet sample was estimated using the 
differences in the transverse and longitudinal shower shapes of multiple photons from meson decays and single photons 
p8|. In the Et range 10 to 50 GeV, 25% ± 25% of the "fake" photons in the background sample were attributed 
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to direct photons. This fraction was subtracted from V^j "7")- The misidentification probabihty was found to 
be V{j "7") - 4 X 10"'' (- 6 X IQ-^) in the CC (EC) in the Et region between 10 and 40 GeV. The measured 
probabihty, before direct photon subtraction, for a jet to mimic a photon is shown in Fig. p7| . 

The total numbers of W + jets background events were estimated to be 1.7 ± 0.9 and 1.3 ± 0.7 for W{ev)^ and 
W{iJLv)"f, respectively, by applying V{j "7") to the observed Et spectrum of jets in the inclusive W{iv) sample. 
The uncertainty on the background estimates was dominated by the uncertainty on V{j "7") due to the direct 
photon subtraction. Bias in the W + jets background estimate due to a possible difference in jet fragmentation 
(e.g. the number of tt^'s in a jet) between jets in the W sample and those in the multi-jet sample was investigated by 
parameterizing ^(j "7") as a function of the EM energy fraction of the jet. No statistically significant difference was 
found between the background estimates with and without the parametrization. The estimated W^+jets background 
also included the background from f+jets, where £ was a jet misidentified as an electron, a cosmic ray muon or a fake 
muon track, since it was derived from the observed inclusive W ^ li' event sample. 

Other backgrounds considered and found to be negligible included those from single photon events where a jet 
was misidentified 
as an electron, and ee + X events where an electron was misidentified as a photon due to tracking inefficiency. 







W{fiu)'y 


Source: 






W -l-jets 


1.7 ±0.9 


1.3 ±0.7 


Zl 


0.11 ±0.02 


2.7 ±0.8 


W{ru)'y 


0.17 ±0.02 


0.4 ±0.1 


Total Background 


2.0 ±0.9 


4.4 ± 1.1 


Data 


11 


12 



TABLE VII. Summary of W{eu)'y and W{fiv)'y data and backgrounds. 



D. Cross section and limits on the coupling parameters 

After subtraction of the estimated backgrounds from the observed number of events, the number of signal events 
was found to be 

JSlWieuh ^ 9_Q+4.2 ^ Q j^Wif.^h ^ 7,Q+f;^ ± 1.1, 

where the first uncertainty is statistical, calculated following the prescription for Poisson processes with background 
given in Ref. |2^, and the second is systematic. 

Using the acceptance for SM couplings of 0.11 ± 0.01 for W{ev)^ and 0.29 ± 0.02 for W{iiv)^ and the efficiencies 
quoted above, the cross section (for photons with Ej, > 10 GeV and ATZe-^ > 0.7) was calculated from a combined 
e + /i sample: 

aiW-f) = 138t^8(stat) ± 21(syst) pb, 

where the systematic uncertainty includes the uncertainties in the e//i/7 efficiencies, the choice of the structure 
functions and the scale at which the structure functions are evaluated, the px distribution of the Wj system, and 
the integrated luminosity calculation. The systematic uncertainties from the sources other than trigger and lepton 



selection efficiencies and geometrical acceptances are listed in Table JVni. The measured cross section agrees with the 
SM prediction of cr^/!^ = 112 ± 10 pb within the uncertainty. Figure |18| shows the data and the SM prediction plus the 
background in the distributions of E^, ATZ^'^, and the three-body "cluster" transverse mass defined by Mt{'j(^; v) — 
[{{m^g + + ExP)^ ±-^t)^ — |Ei^ + ±!pTp]^- Final state radiation events and background events composed 
most of the expected signal with MT{'^t',v) < Mw- Of the 23 observed events, 11 events had MT{'^t',v) < Mw- 
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Uncertainty 


Luminosity 


5.4% 


Structure Function choice 


6.0% 




3.9% 


Conversion Probability 


5.0% 


Random track overlap 


1.0% 


Photon selection efficiency 


7.0% 


Total 


12.5% 



TABLE Vlll. The values of systematic uncertainties in the cross section and coupling limit measurements, other than 
those of trigger, lepton selection and acceptance. 



To set limits on the anomalous coupling parameters, a binned maximum likelihood fit was performed on the E'^ 
spectrum for each of the W{ev)^ and W{^v)^ samples, by calculating the probability for the sum of the Monte 
Carlo prediction and the background to fluctuate to the observed number of events (see Appendix 2 for more detail). 
The uncertainties in background estimate, efhciencies, acceptance and integrated luminosity were convoluted in the 
likelihood function as Gaussian distributions. A dipole form factor with a scale A = 1.5 TeV was assumed for the 
anomalous couplings in the Monte Carlo event generation. The Monte Carlo events were generated at 11 x 11 grid 
points of the CP-conserving anomalous coupling parameters, Ak^ and A^, assuming that the CP-violating anomalous 
coupling parameters and A-y are zero. The limit contours for Ak^ and A^ are shown in Fig. The numerical 
values of the limits at the 95% confidence level (CL) were 

-1.6 < Ak^ < 1.8 (A.^ = 0), - 0.6 < A.^ < 0.6 (Ak^ = 0) 

for s = (i.e. the static limit). The J7(l)£;M-only coupling of the W boson to a photon, which leads to K-y = (Ak^ = 
— 1) and A7 = 0, and thereby, = e/2mw and = [Q, was excluded at the 80% CL, while the zero magnetic 
moment {^w — 0) was excluded at more than the 95% CL. Similarly, limits on CP-violating coupling parameters 
were obtained as —1.7 < < 1.7 (A-y = 0) and —0.6 < A^ < 0.6 {k^ = 0) at the 95% CL. The form factor scale 
dependence of the results was studied. It was found that the limits were insensitive to the values of the form factor 
scale for A > 200 CeV and were well within the constraints imposed by S-matrix unitarity |3^] for A = 1.5 TeV. A 
simultaneous fit to and the A7?.^7 spectra was performed. It was found that the results were within 3% of those 
obtained from a fit to the E^ spectrum only. 



VI. WW DILEPTONS 

In this section the results of a search for pp — > WW + X ^ Uvv' + X, where the leptons included muons and 
electrons, are presented. The signal and background were estimated and an upper limit was set for the cross section 
of the SM process. Anomalous WWZ and WW^ couplings would have enhanced the expected WW cross section 
by upsetting the cancellation |Q between the production diagrams and the trilinear diagram as seen in Fig. IJ, 
which shows the cross section vs anomalous couplings for A = 1000 GeV. The detection efficiency also increases with 
anomalous couplings because of the higher average Ex of the W bosons (see Fig. ||), resulting in a higher average Ej^ 
for leptons and in more central events. This expected increase in the cross section and efficiency was exploited to set 
limits on the anomalous coupling parameters, A and An. 

The expected signature for W boson pair production with subsequent decay to dileptons was two high-p^ isolated 
leptons in association with large I^t- The major sources of background were the following: events with a W + jet(s) 
where a jet was misidentified as a lepton; W + ^ events where a photon was misidentified as an electron; QCD multi-jet 
events where two jets were misidentified as leptons; Z ^ U, Z ^ tt ^ ll'vvvv events; and tt ^ + X events. The 
event selection requirements were designed to reduce these backgrounds while retaining high detection efficiency for 
signal events. The selection requirements were slightly different for the ee, e^, and /x/i channels because the electrons 
had a better pj- resolution but a larger background contamination than muons. In what follows, the analyses of 
individual channels and our limits on the cross section for W boson pair production as well as on the anomalous 
gauge boson trilinear couplings are presented. 
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A. The ee channel 



The WW eevv candidate events were selected from the data sample recorded using the ELE-2-HIGH trigger 
which required two EM clusters with Et > 7 GeV at level 1 and two isolated EM clusters with Et > 10 GeV at 
level 2 (see Table ^. Candidate events containing two electrons that passed the "Tight I" requirements were selected. 
The "Tight I" requirements discussed in Section ||, and detailed in Table III, provided the largest rejection of fake 
electrons. The following event selection requirements were then imposed. Both electrons were required to have a large 
transverse energy {Et > 20 GeV); at this stage the remaining sample of 605 events was comprised primarily of Z 
bosons. The I^t of the event was then required to be > 20 GeV. These first two selection criteria strongly reduced the 
background due to QCD fakes. The dielectron invariant mass was required to be outside of the Z boson mass window 
(between 77 and 105 GeV/c^). The I^t and dielectron invariant mass selections had very strong rejection (> 100) 
of Z ee decays. The background from Z ^ tt eevv which was not eliminated by the electron pT thresholds 
was further reduced by requiring that the Ij^T not be coUinear with the direction of the lower energy electron; it was 
required that 20° < A(/)(p|n, _^t) < 160° for the lower energy electron if Ifx < 50 GeV. Releasing this requirement 
for events with large Ifx increased the acceptance for W boson pairs in a region where the Z boson background was 
very small. This can be seen in Fig. ^ which shows A(/)(p|,, I^t)^s I^t distributions for W boson pairs, Z ^ ee and 
Z ^ rr — > eevv. Finally, the sum of the Et of the recoiling hadrons (E'^*^'*), defined as —{E!^ + El^ + I^t) was 
required to be less than 40 GeV in magnitude. The background from ti production was effectively eliminated by this 
requirement. Figure |l| shows a Monte Carlo (PYTHIA plus D0GEANT) simulation of E^""^ for - 20 fb^^ of SM 
WW and tt events. For WW events, non-zero values of E^^'^ were due to gluon radiation and detector resolution. 
For tt events, the most significant contribution was from h quark jets from t quark decays. This selection reduced 
the background from ti production by a factor of more than four for a t quark mass of 170 GeV/c^. The efficiency of 
this selection criterion for SM W boson pair production events was 0.951qq4 and decreased slightly with increasing 
W boson pair invariant mass. The systematic uncertainty in the efficiency of this last selection criteria, included in 
the uncertainties presented, was estimated from the difference between E^^'^ for Z boson data and Monte Carlo ( 
PYTHIA plus D0GEANT) distributions. 

Table IX shows the numbers of events remaining after each selection cut. One event survived all the selection 
criteria. 



Event selection criteria Number of events surviving 
Et > 20 GeV 605 
^t > 20 GeV 5 
Mee < 77GeV/c^ or M^^ > 105 GeV/c^ 3 
A0(p5?,^t) cut 1 
0had^ < 40 GeV 1 



TABLE IX. The numbers of events after each selection cut for the WW ee analysis. 



Fiducial Region 


Efficiency 


cc - cc 


0.526 ±0.041 


CC-EC 


0.368 ± 0.044 


EC - EC 


0.257 ±0.058 


TABLE X. The combined trig 


ger and electron selection efficiency for individual fiducial regions in the WW ^ ee analysis 


for SM W pair production. 





The integrated luminosity of the data sample was 14.3 ± 0.8 pb^^. The trigger efficiency was calculated with 
the TRIGSIM simulation package to be etrig — 0.989 ±0.002. The measured electron selection efficiency was used 
to estimate the detection efficiency for SM WW eevv events. The geometrical acceptance was obtained from 
a PYTHIA and D0GEANT Monte Carlo simulation. These efficiencies for individual fiducial regions are listed in 
Table The overall detection efficiency for SM W pair events was estimated to be Eee — 0.094 ±0.008. The expected 
number of events N^^^ was N^^^ — 0.149 ± 0.013(stat) ± 0.019(sys), using the next-to-leading order cross section 
||3|, and branching fraction Br(I^ ev){= 0.108 ± 0.004) g]. 

The backgrounds from M^7, Drell-Yan dilepton, Z ^ tt ^ eevvvv and ti processes were estimated using the 
PYTHIA and ISAJET Monte Carlo event generators followed by the D0GEANT detector simulation. The ti cross 
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section estimates were from the calculations of Laenen et al. [Q. The tt background was averaged for Mtop = 160, 
170, and 180 GeV/c^. The production and decay of Z bosons was modeled using the double differential (in rapidity 
and pt) cross section calculated at next-to-leading order Q, and a fast detector simulation of the type discussed in 
Section [II A. The line shape of the Z boson was taken to be a relativistic Breit-Wigner function. The kinematic 
distributions were compared with the Z boson data sample and found to be consistent. 

The backgrounds from W^-l-jet(s) with a jet misidentified as an electron and multijet events with two jets misidentified 
as electrons were called "fake" background. The size of this background was estimated with the following method. 
Two sub-samples of data were derived from the full data set. One was similar to the signal sample and contained two 
"Tight I" electrons each with Et > 20 GeV. The other was a sample of events with at least one bad electron which had 
an H- matrix > 200 and Isolation fiso > 0.15 (the fake sample). A normalization factor (Ffake) of this fake sample 
relative to the signal sample was calculated using the number of events with I^t < 15 GeV, which contained solely fake 
electrons, in both the samples. All the event selection cuts were applied to the fake sample and the number of remaining 
events [Nfake] was counted. The fake background {Nj^J was then computed from the products of Ffake and Nfake- 
The result was N^^^ = 0.152 ± 0.012(stat) ± 0.076(sys). The total number of background events was estimated to be 
Ng^ = 0.222 ± 0.020(stat) ± 0.080(sys). Table g contains a summary of the expected background in the ee channel. 



Background 


ee 


e/i 




Z ^ ee or fifi 


0.02 ±0.01 




0.068 ± 0.026 


Z ^ TT 


< 10"=* 


0.11 ±0.05 


< 10"^ 


Drell-Yan dileptons 


< 10"^ 




< 10"^ 


Wy 


0.02 ±0.01 


0.04 ±0.03 




QCD (iV/„t) 


0.15 ±0.08 


0.07 ±0.07 


< 10"^ 


tt 


0.03 ±0.01 


0.04 ±0.02 


0.009 ± 0.003 


Total 


0.22 ±0.08 


0.26 ±0.10 


0.077 ±0.026 



TABLE XL Summary of the expected number of background events to WW — > eei/u, WW ej-ii^u and WW n^vu. The 
uncertainties include both statistical and systematic contributions. 



B. The e/i channel 

The WW efivv candidate events were selected from the data sample acquired with the MU-ELE trigger 
which required at least one EM tower with Et > 7 GeV and one muon with I77I < 1.7 at level 1 and one EM 
cluster with Et > 7 GeV and one muon with pT > 5 GeV/c at level 2 (see Table ^. The electron in a can- 
didate event was required to pass "Tight I" criteria, providing the strictest rejection against fake electrons, and 
the muon to pass the "Loose I" criteria of Table |l|. The following event selection requirements were imposed. 
Both the electron and the muon were required to have a large transverse energy (momentum) E^ > 20 GeV 
and pt^ > 15 GeV/c. Both I^t and of the event were required to be > 20 GeV. These first two require- 

ments provided large rejection of the background from multijet events. In order to reduce the backgrounds from 
Z ^ TT ^ efivvvv, the I^t was required not to be coUincar to the muon: 20° < A(j){pt^, I^t) < 160° if I^t < 50 
GeV. Figure |22| shows the A0(p^,^t) vs I^t distributions. Finally, the recoil hadronic Et {Et^'^)^ defined as 
— {E^ + EI^+JPt) was required to be less than 40 GeV in magnitude to reduce the background from ti produc- 



tion. Table XI] shows the number of events remaining after each selection cut. One event survived all the require- 



ments but the last; this event is a candidate for tt production and has been discussed extensively elsewhere [p5| 



Event selection cut 


Number of events 


> 20 GeV 


9 


pt^ > 15 GeV 


6 


> 20 GeV 


1 


A0(j3^,^t) cut 


1 








TABLE XII. The numbers of events remaining after each selection criteria for the WW —> ev^v analysis. 
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The integrated luminosity of the data sample was 13.9 ± 0.8 pb~^. The trigger efficiency was largely deter- 
mined by the trigger efficiency for the muons and was estimated using data as was discussed in Section ||. The 
detection efficiency for the SM WW e^vv events, including the muon selection efficiency, the geometrical ac- 
ceptance, and the event selection efficiency, was estimated using the PYTHIA and D0GEANT Monte Carlo sim- 
ulation. The muon selection efficiency was implemented in the D0GEANT Monte Carlo program by introduc- 
ing the measured hit efficiencies and resolutions of the muon chamber modules. The measured electron selec- 
tion efficiency was implemented as a multiplicative factor after the detector simulation. The detection efficien- 
cies, including lepton identification efficiencies, of individual fiducial regions are listed in Table XIII. The un- 
certainty on the efficiency for the regions involving EF muons was dominated by the statistics of the GEANT 
Monte Carlo simulation. The overall detection efficiency of the SM W pair events was estimated to be e^^ — 
0.092 ± 0.010. The expected number of events was estimated to be N^^' = 0.283 ± 0.031 (stat) ± 0.037(sys). 



Fiducial Region 


Efficiency 


CC - CF 


0.43 ±0.08 


CC - EF 


0.21 ±0.14 


EC - CF 


0.30 ± 0.09 


EC - EF 


0.15 ±0.15 



TABLE XIII. The combined efficiencies of trigger, electron and muon selection, and kinematic event selection for individual 
fiducial regions in the WW — > envu analysis. The fiducial regions are those of the charged leptons. 



The backgrounds from W"j, Z ^ tt ^ eyuvvvv and tt were estimated using the PYTHIA and ISAJET Monte Carlo 
event generators followed by the D0GEANT detector simulation. The background due to a jet misidentified as an 
electron was estimated by a different method from the ee channel, since the accuracy of the estimate was limited 
by statistics when that method was applied to the e/x channel. Instead, the inclusive W iiv data were used to 
estimate this background. Each jet in an event was treated as an electron and the event selection requirements were 
applied. The events that survived the criteria were weighted by the probability of a jet being misidentified as an 
electron. The misidentification probabilities were measured from data to be P'^'^(jet — > e) = (0.9 ± 0.4) x 10"'' for 
CC and P^'-^(jet ^ e) = (4.0 ± 1.0) x lO"** for EC. The total "fake" background from this source was calculated to be 
N^g = 0.074 ± 0.016(stat) ± 0.074(sys). The background due to a jet faking a muon was estimated to be negligibly 
small. The total number of background events was estimated to be N^^ = 0.264 ± 0.052(stat) ± 0.084(sys). The 



backgrounds to WW efivv are summarized in Table XI 



C. The nfj. channel 



The WW jjijjiuu candidate events were selected from the data sample recorded with the MU-MAX trigger of 
Table |. This trigger required at least one muon in |?7| < 1.7 at levels 1 and 1.5 with pT > 7 GeV/c (threshold 
determined by the hardware) and one muon with px > 15 GeV/c at level 2. Candidate events with two muons that 
passed "Tight I" requirements were selected. The following event selection requirements were imposed. Both muons 
were required to have large transverse momentum: > 20 GeV/c and pf^ > 15 GeV/c. To reduce the background 
from Z ^fi decays, the was required to be > 30 GeV, where was defined as the projection of the I^t vector 
onto the bisector of the opening angle of the two muons in the transverse plane. By selecting this component of 
the I^T, it was ensured that the I^t was least sensitive to mismeasurements of the muon momentum. This selection 
requirement was also less sensitive to the momentum resolution of the muons than was a dimuon invariant mass 
cut. The I^T was required not to be coUinear to the higher momentum muon: A0(p^\^t) < 170°. This reduced 
the background from Z tt —f p.^vvvv. Figure [2^ shows A(j){pt^ , I^t) vs I^t distributions. The recoil hadronic Et 
{E}^'"^), defined as -{E^'^ + Elf + $t) was required to be less than 40 GeV in magnitude, rejecting tt as in the other 
two channels. Even though the pr resolution of the muons was worse than that of the electrons, the resolution of this 
variable was the same in all three channels since the mismeasurement of the leptons cancels when taken in a vector sum 
with the IpT- Table XIV shows the numbers of events after each selection cut. No event survived all the selection cuts. 
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Event selection cut 


Number of events 


pt^' > 20 GeV/c 


102 


ptj? > 15 GeV 


88 


^t" > 30 GeV 





A0(p^\^t) < 170° 











TABLE XIV. The number of events remaining after each selection criteria for the WW — > ^^vv analysis. 



Fiducial Region 


Efficiency 


CF - CF 


0.023 ± 0.006 


CF - EF 


0.009 ± 0.002 


EF-EF 


0.0010 ± 0.0006 



TABLE XV. The combined efficiencies of trigger and muon selection for individual fiducial regions in the WW fifivu 
analysis. CF-CF is, for instance, the case that both muons were in the central region. 



The integrated luminosity of the data sample was 12.2 ± 0.7 pb^^. The trigger efficiency was measured using data. 
The detection efficiency for the SM WW — > iifiiyjy events that included the muon selection efficiency, the geometrical 
acceptance and the event selection efficiency was estimated using the PYTHIA and D0GEANT Monte Carlo simulation. 
The muon selection efficiency was implemented in the D0GEANT Monte Carlo program by introducing the measured 
hit efficiencies and resolutions of the muon chamber modules as in the e/i channel. The efficiencies of individual 
fiducial regions are listed in Table XV. The overall detection efficiency for the SM W pair events was estimated to 
be Cf^f, = 0.033 ± 0.006. The expected number of events was estimated to be N^f^'^ = 0.045 ± 0.004(stat) ± 0.006(sys). 

The backgrounds from Drell-Yan dilepton, Z ^ tt ^ jifivvvv and tt processes were estimated using the PYTHIA 
and ISAJET Monte Carlo event generators followed by the D0GEANT detector simulation. The background from 
Z ^/z was estimated using the same fast simulation program as in the ee channel. The fake background due 
to a jet faking a muon was negligibly small. The total number of background events was estimated to be — 
0.077 ± 0.023(stat) ± 0.012(sys). 



D. Limit on the cross section for W boson pair production 



The results from the analyses of the ee, e/i and /x/i channels are summarized in Table XVI. For the three channels 
combined, the expected number of events for SM W boson pair production, based on a cross section of 9.5 ±1.0 
pb , was 

0.47±0.07. In approximately 14 pb~^ of data, one event was found with an expected background of 0.56±0.13 events. 



Channel 


ee 




MM 


Total 


Efficiency 


0.094 ±0.008 


0.092 ±0.010 


0.033 ± 0.006 






0.15 ±0.01 ±0.02 


0.28 ± 0.03 ± 0.04 


0.045 ±0.004 ±0.006 


0.47 ±0.03 ±0.06 




0.22 ±0.02 ±0.06 


0.26 ±0.05 ±0.08 


0.077 ± 0.023 ±0.012 


0.56 ±0.06 ±0.10 


jt^observed 


1 








1 



TABLE XVI. The summary of WW dileptons analyses including the efficiency, number of SM events expected, expected 
backgrounds, and number of candidates observed. 



The 95% confidence level upper limit on the W boson pair production cross section was estimated based on one 
observed event, taking into account the expected background of 0.56 ± 0.13 events. Poisson-distributed numbers of 
events were convoluted with Gaussian uncertainties on the detection efficiencies, background and luminosity. For SM 
W boson pair production, the upper limit for the cross section was 87 pb at the 95% confidence level. 
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E. Limits on the trilinear gauge boson couplings 

The limit on the W boson pair production cross section can be translated into limits on the anomalous gauge boson 
couplings. The Monte Carlo program of Ref. followed by a fast detector simulation was used to estimate Q the 
detection efficiency for W boson pair production as a function of the coupling parameters A and Ak. It was assumed 
that the W boson couplings to the photon and to the Z boson were equal: A = A-y = and Ak = Ak^ = Akz ■ The 
form factor scale A = 900 GeV was chosen. This was the highest value of A that produced anomalous coupling limits 
within the corresponding unitarity bound in this analysis. For smaller values of A, the anomalous coupling limits are 
looser (see Appendix 1). The MRSD— ' parton distribution functions were used in the event generation. 

The number of W boson pair events expected at each point in a grid of A and k, including SM production, was 
fitted to the following equation which reflected the general Lagrangian form of gauge boson self-interactions: 

N{Ak,X)= ai + a2AK + a3(AK)2 
+ a^X + a5(A)^ + flgAAK, 

where the Ui were parameters determined from the fit. A 95% CL limit contour on the coupling parameters Ak and 
A was formed by intersecting the parabolic surface of expected number of events with the plane of the 95% CL upper 
limit on the observed number of events (with the background subtracted) , fluctuated by Gaussian uncertainties on the 
detection efficiencies, backgrounds, and luminosity and by the Poisson uncertainty on the statistics of the observation. 
The 95% CL upper limits on the coupling parameters are shown in Fig. ^ (solid line). Also shown in Fig. ^ (dotted 
line) is the contour of the unitarity constraint on the coupling limits for the form factor scale A = 900 GeV. This 
value of A was chosen so that the observed coupling limits lie within this ellipse. The limits on the CP-conserving 
anomalous coupling parameters were —2.6 < Ak < 2.8 (A = 0) and —2.1 < A < 2.1 (Ak — 0). The limits for the 
CP-violating parameters, k and A, were similar. The limits on A and Ak exhibited almost no correlation, in contrast 
to limits from the W'y analyses presented in Section ^ and in Refs. |ll| , [l^ . 



VII. SEARCH FOR ANOMALOUS WW AND WZ PRODUCTION IN THE ei^ JET JET CHANNEL 



In this section, a search for anomalous WW and WZ production is presented. The method was to identify WW 
and WZ candidates where one W boson decayed to an electron and a neutrino and the Z boson or other W boson 
decayed to two jets. The expected cross section times branching fraction for the SM WW and WZ processes (^ 1.6 pb) 
was much smaller than that expected from the W boson plus dijet background (~ 76 pb) |37|, which had similar 
characteristics. Rather than isolating the SM signal, limits were set on the anomalous couplings by comparing the 
characteristics of the events with those expected from non-SM couplings. Figure ^ shows the pT of the W bosons 
for SM WW production and for WW production with anomalous trilinear couplings. Anomalous couplings lead to a 
dramatic increase in the cross section at high pt{W). To exploit this, the pr of W bosons in the candidate events was 
measured, the contribution of the backgrounds to that spectrum was estimated, and the data were compared with 
the sum of the background plus the expectations for the signal for various anomalous couplings. 



A. Event Selection and Efficiency 



The WW {W Z) ei/jj candidates were selected by searching the data which passed the ELE-HIGH trigger (see 
Table |) for events with a high Et electron accompanied by significant I^t and at least two jets consistent with 
W jj or Z — > jj. Events with electrons which satisfied the "Tight IV" criteria within \ri\ < 2.5 and with Et > 25 
GeV were chosen. The I^t was required to be greater than 25 GeV, and at least two jets were demanded, each with 
Et > 20 GeV and \ri\ < 2.5. A small jet cone size, TZ < 0.3, was used to ensure that the two jets from the W or 
Z decay, close together for high-p^ W or Z bosons, were resolved into distinct jets. After the jets were identified, a 
cleanup algorithm was applied to remove events with "fake" jets due to noisy cells or badly mismeasured jets, which 
occured primarily in the intercryostat region. Figure ^ shows the transverse mass of the electron and I^t, M^'^, 
for the candidates which survived the preceding selection criteria. MjVwas required to be greater than 40 GeV/c^. 
The dijet invariant mass distribution of these events is shown in Fig. Eq. In case there were more then two jets with 
Et > 20 GeV in the fiducial region, the combination yielding the largest invariant mass was taken to be the dijet 
mass of the candidate W oi Z boson. Requiring the dijet invariant mass to be 50 < mjj < 110 GeV/c^ yielded 84 
candidate events. 
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The efficiency for identifying two separated jets depended on the pr of the W boson. For pt{W) < 125 GeV/c, 
the efficiency was dominated by the jet Et threshold. For pt{W) > 350 GeV, the efficiency was dominated by the 
probability for the two jets to merge into one in the reconstruction process (hence the use of the small cone size). Using 
the ISAJET and PYTHIA event generators, followed by the detailed detector simulation, D0GEANT, and the Shower 
Library described in Section III A, the efficiency for reconstructing W jj was estimated as a function of pt{W), 
including the jet-finding efficiency and the efficiency for the dijct mass requirement. The Z ~f jj reconstruction 
efficiency was obtained in a similar manner. From the Monte Carlo it was determined that the use of the two highest 
Et jets to form the dijet mass was the correct assignment 90% of the time. Figure ^ shows the efficiency for the dijet 
reconstruction of jj as a function of pt{W) for events generated with ISAJET and PYTHIA. For the efficiency, 

the results from the ISAJET simulation were used because they were smaller than the efficiencies determined from 
PYTHIA; the difference (9%) was included in the systematic uncertainty for the efficiency. 

The overall efficiency was calculated for SM and anomalous couplings using the fast detector simulation de- 
scribed in Section II] along with the WW {WZ) generator |jl^. The reconstruction efficiencies for W and Z 
boson decays to dijets were incorporated as look-up tables. The pT distribution of the WW and WZ sys- 
tems was included in the simulation by using the observed pt{Z) spectrum from the inclusive Z — > ee sam- 
ple. The uncertainty in the absolute jet energy scale and in the jet corrections, 10%, was included in the 
systematic uncertainty by recalculating the results shifting the jet energies within their uncertainty. Other 
sources of uncertainty included: 6% for the uncertainty in the resolution of the I^t and 4% for the uncer- 



tainty in the electron identification efficiency. All of the uncertainties in the efficiency are listed in Table XVII 



Source 


Uncertainty (%) 


Statistical 


1 


Electron Efficiency 


4 


l^T Smearing 


6 


Jet Energy Scale 


6 


Jet Reconstruction Efficiency 


9 


Total 


13 



TABLE XVII. Summary of systematic uncertainties for the WW/WZ — > eujj analysis. 



The total efficiency for the detection of SM WW and WZ events was estimated to be 0.15 ± 0.02 and 0.16 ± 0.02, 
respectively. Therefore the total number of expected (SM) signal was 3.2 ± 0.6 events where 2.8 ± 0.64 events were 
WW and 0.4 ± 0.1 events were WZ, including the uncertainties in the efficiency and luminosity. 



B. Background Estimate 



The background included contributions from the following: W+ > 2 jets; ti production with subsequent decay to 
W^W^bb, where the top mass was assumed to be 180 GeV/c^; WW{WZ) production with W tv followed by 
r — » evD] ZX — > eeV, where one electron was mismeasured or not identified; and multijet events, where one or more 
jets was misidentified as an electron and there was significant I^t due to mismeasurement or the presence of neutrinos. 

The multijet background was estimated following the same procedure used in the WW — > eevv channel. The 
background sample was comprised of events which contained a jet with an EM fraction greater than 0.9 within the 
electron fiducial region and a matching track. However, these electron candidates satisfied at least one of the following 
three electron "anti-identification" criteria: fjso > 0.15, H-matrix > 250, or track match significance, TMS > 10. 
The number of events in the region < I^t < 15 GeV/c was used to normalize the fake sample to the signal sample. 
This was done after all selections except for the dijet mass cut. Then the dijet mass selection criterion was applied to 
the fake sample to determine the number of background events. The possible signal contamination of the fake sample 
was included as a 7% systematic uncertainty in the normalization of the fake background. A systematic uncertainty 
of 3.4% comes from the variation in the fake event normalization when the upper end of the normalization region was 
varied in the range 12 to 18 GeV. An uncertainty of 4% arose from variation of the I^t threshold of the signal over 
the range 22 to 28 GeV. Figure ^ shows the I^t for the QCD background and for the signal candidates before the 
dijet mass selection. 

The backgrounds from tt, WW, WZ rvjj and ZX — > eeX were estimated from D0GEANT simulation of 
PYTHIA and ISAJET events. The background from W+ > 2 jets came from VECBOS Q generated events car- 
ried through a hybrid detector simulation which combined the D0GEANT detector simulation with the parton- 
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based jet Shower Library. The normahzation of the W+ > 2 jets background, made before the dijet mass se- 
lection, was determined from the comparison of the number of candidate events outside the dijet mass window 
with that expected from the VECBOS Monte Carlo, after subtracting the multijet backgrounds and expected 
SM signal. The systematic uncertainty in the W+ > 2 jets background included contributions from uncertainty 
in the fake normalization amounting to 7%; variation when the dijet mass window was increased in width to 
40 < rrijj < 120 GeV/c^, amounting to 10%; and variation in the background when the Monte Carlo jet en- 
ergy scale was increased by 10%, amounting to 11%. The cross section for the resulting W+ > 2 jets back- 



ground agreed within 1.5% of the VECBOS expectation. The background estimate is summarized in Table XVIII 



Background Source 


N 

Before Dijet Mass Cut 


amber of Expected Events 

After Dijet Mass Cut 


W+ > 2 jets 
tt 

WW, WZ ^ Tujj 
ZX eeX 
multijets 


125.4 ± 25.9 
3.42 ± 0.47 
0.24 ±0.02 

"•'-"-'-0.00 

30.0 ±4.5 


62.2 ± 13.0 
0.87 ±0.12 
0.22 ±0.02 

12.2 ± 2.6 


Total Background 


159 ± 26 


75.8 ± 13.3 


SM WW + WZ Prediction 


3.4 ±0.6 


3.2 ±0.6 


Data 


166 


84 



TABLE XVIII. Summary of WW(WZ) evjj backgrounds and data. 



The distributions in pxieiy) of the final event sample, and for the W+ > 2 jets background, the total background, 
the SM WW and WZ Monte Carlo, and the WW and WZ Monte Carlo for a non-SM value of the couplings 



{Akz — — 2, \z = X-y = 1.5) are shown in Fig. 29, The px spectrum was consistent with that expected from the 
background. The highest pt event, important in setting anomalous coupling limits, had pT{ev) ~ 186 ± 20 GeV/c. 
There were no other candidates with pT{ev) > 100 GeV/c. 

Using the detection efficiencies for SM WW and WZ production and the background-subtracted signal, and as- 
suming the SM ratio of cross sections for WW and WZ production, an upper limit at the 95% confidence level (CL) 
on the cross section a{pp W~^W~X) of 183 pb was determined. 



C. Determination of Limits on Anomalous Couplings 

The absence of an excess of events with high pt{W) excluded large deviations of the trilinear couplings from the 
SM values. The px spectrum expected at each element in a 225 point grid in A and An space, centered around and 
including the SM values, was obtained using Monte Carlo simulation. The assumptions on the A scale and on the 
relation between the WW'y and WW Z couplings affected the pt spectra. Unequal width bins were used in order to 
evenly distribute the events, particularly at the ends of the spectrum. An analytic form of the prediction of the number 
of events in each pt bin was obtained with a quadratic function of the coupling parameters, similar to that used in 
the dilepton analysis, fit to the number of events for each pair of anomalous couplings. The difference between the 
estimated number of events and the fit was calculated for each pair of anomalous couplings for a particular pT bin and 
found to be less than 10%. This value was included as a systematic uncertainty in the fitting procedure. To set limits 
on the anomalous couplings, a binned likelihood fit was performed on the pt(W^) spectrum of the expected signal 
plus background for pt(W^) > 25 GeV/c. In each pT bin, the probability was calculated for the predicted number of 
events to fluctuate to the observed number of events. The uncertainties in the efficiency, background estimates, and 
total luminosity were convoluted in the likelihood function using Gaussian distributions. This likelihood fit procedure 
is described in detail in Appendix 2. 

Limits were obtained on the coupling parameters under four sets of assumptions on the relations among the 
coupling parameters. For all four assumptions, the most likely point in the A — Ak grid was the SM point. For 
the assumption Ak = Are-y = Ak^, and A = A-^ = A^ with A = 1500 GeV, the contours for the 95% CL limit 
on A and Ak, with A = 1500 GeV, are shown in Fig. |^(a). The 95% CL hmits were -0.9 < Ak < 1.1 (A = 0) 
and —0.6 < A < 0.7 (Ak = 0). As in the WW dileptons analysis, the limits on A and Ak exhibited almost 
no correlation. Under the HISZ relations |^^, which parameterize the WWZ couplings in terms of the WWj 
couplings: Akz — 0.5 Ak.^ (1 — tan^ 6'tu), Agz — 0.5 Ak^/ cos^ 9^, Xz — X^, the 95% CL coupling limit contours with 
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A = 1500 GeV are shown in Fig. gO|(b). The Hmits were -1.0 < Ak^, < 1.3 (A^ = 0) and -0.6 < A..^ < 0.7 (Ak^ = 0). 
Under the assumption that the WW^ couphngs have the SM value, the 95% CL upper hmit contour, in Xz and Akz, 
is shown in Fig. |o|(c). The 95% CL hmits were -1.1 < Akz < 1.3 {Xz = 0) and -0.7 < Xz < 0.7 {Akz = 0). 
Under the assumption that the WWZ couphngs have the SM value, the 95% CL upper limit contour, in A-^ and Akj, 
is shown in Fig. |o|(d). Here the A scale was 1000 CeV. The 95% CL Hmits were -2.8 < Ak^ < 3.3 (A.^ = 0) and 
—2.5 < A^ < 2.6 {Ak^ = 0). The limits from 5-matrix unitarity are also shown in Fig. ^(a)-(d) for each assumption. 
The unitarity limits were ellipses for (a) and (b) due to the form of Equation 11.1, shown in Appendix 1. However, 
for (c) and (d), the intersections of the Wj and WW/WZ unitarity contours are shown in the Figure. 

Because this analysis accounted for the background in fitting the spectrum for pt{W) > 25 GeV/c, it was sensitive 
to anomalous couplings at both large and small s. All of the results of the fits were insensitive to the priW) threshold 
when varied between 25 and 130 GeV/c. In contrast, the analysis in Ref. |l^, which required pt{W) > 130 GeV/c, 
loses sensitivity at small s; deviations from the SM restricted to s < 500 GeV could have been missed pO|. 



VIII. COMBINED W-y AND WW/WZ ANOMALOUS COUPLING RESULTS 

The WW dileptons counting experiment and WW/WZ pT spectrum analysis are sensitive to the same WW^ 
couplings as the W^ photon spectrum analysis. The three analyses can be combined to form tighter limits on 
anomalous couplings. In this section, the procedure and results of the combined fit are discussed. 

The likelihood method used in the WW/WZ pT analysis and W^ photon spectrum analysis was used in the 
combined analysis. The joint log-likelihood was the sum of the log of the probabilities, as discussed in Appendix 2. 
The likelihood was formed from the Monte Carlo WW/WZ pr spectrum and Monte Carlo Wj photon spectrum, 
expected background, and observed number of events in each channel with identical binning as was used in the 
separate analyses. The expected number of WW dilepton events was recalculated for A = 1500 GeV (equivalent 
to use of a single bin for all pt{W)); while the A and Ak limits would have violated unitarity for this value of A, 
the combined limit does not. Common systematic uncertainties, including lepton identification efficiency (4% for 
all channels with an electron in the final state and 12% for all channels with a muon in the final state), integrated 
luminosity (5.4%), and choice of parton distribution function (9.1%), were treated as discussed in Appendix 2. The 
limits are insensitive to a change in the size of the common systematic uncertainty by as much as a factor of two. The 
statistical uncertainties of the data dominate the uncertainty in the analysis. 

The following results were obtained. For the assumption that the WW^ couplings are equal to the WWZ couplings 
and with A = 1500 GeV, the 95% CL limits were -0.71 < Ak < 0.89 (A = 0) and -0.44 < A < 0.44 {An = 0). 



Figure 31 shows the 95% CL limit contour for A and Ak along with the unitarity contour. 



IX. ^7 PRODUCTION 

A measurement of the ZZj and .2^77 couplings using pp — > ££"/ + X {£ ^ e, ^, i^) events is discussed in this section. 

The signature for Zj events was two high-pT leptons (e+e^, fJ-^fJ-^ or vv), and a photon. The leptons would not 
necessarily have combined to give the Z boson mass. In initial state radiation and anomalous coupling events, of the 
type shown in Figs. 0(a)-(c), the dilepton invariant mass for the electron and muon decay channels would be at the 
Z boson mass. However, for events with bremsstrahlung radiation from a charged lepton, as shown in Fig. H, the two 
leptons would have a pair mass below that of the Z boson. Furthermore, photons radiated from the leptons would 
have tended to be close to the leptons. The neutrino decay channels had several important differences. Besides having 
a higher branching fraction than the electron and muon decay channels (20.0% for three generations of neutrinos vs. 
3.37% for ee or /i/i), the Z decays are inferred with high efficiency in the detector through the I/!t measurement. 

The radiative diagrams do not contribute to the neutrino decay channel. Thus, the cross-section changes more quickly 
with anomalous couplings than the cross section for the electron and muon channels. The signature for these events 
was a photon recoiling against the 1/!t of the undetected neutrinos. The main disadvantage of the neutrino channel 
was that the backgrounds were larger than in the other channels. 



A. The ee channel 



The 667 sample was selected from events which satisfied the ELE-2-MAX trigger described in Table ^ The data 
set corresponded to an integrated luminosity of 14.3 ± 0.8 pb^^. From this sample, candidate events were required 
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to have two electrons in the fiducial region with Et > 17 GeV. At least one electron had to satisfy the "Tight 11" 
requirements (see Table HI) while the other satisfied the "Loose I" requirements. This combination of tight and 
loose electron selection was possible because the backgrounds from fake electrons were small, relative to the expected 
signal, when the photon was required to pass the "Loose" requirements (see Table IV) within the fiducial region. After 
trigger, fiducial region, and particle selection criteria were applied, 10 events with two electrons with Et > 17 GeV 
and a photon with Et > 5 GeV survived. Final selected events were required to have E^ > 25 GeV and a photon 
separated from each electron by ATZ^j > 0.7 with E'^ > 10 GeV. These last two requirements reduced the contribution 
of radiative events. Four events survived in the final sample. Table XIX indicates the number of events surviving the 
last few selection criteria. For details on the characteristics of individual events and for event displays see Ref. . 



Selection Criteria 


No. of Surviving Events 


Starting Sample 


77 


Fiducial and Particle ID 


15 


Trigger Criteria 


10 


Et'' > 25 GeV 


10 


A7^e^. > 0.7 


8 


> 10 GeV 


4 



TABLE XIX. Number of eey candidates which passed the selection criteria. 



The trigger efficiency for SM production was estimated using the Z ^ ee event sample. It was found to be 
0.98 ± 0.01. The acceptance for SM Z7 production and for production via anomalous ZZj and Zjj couplings was 



estimated using the event generator of Ref. Il8l combined with the fast detector simulation discussed in Section III 



MRSD— ' structure functions |g7|| were used in the event generation and the cross section was scaled by a k- factor 
of 1.34. The geometric acceptance for SM production was 53%. Averaged over Et for SM production, the photon 
identification efficiency was also 0.53 ± 0.05. With the particle identification criteria, the kinematic, and the fiducial 
requirements on the electrons and photons described above, the selection efficiency for SM Z7 production was 0.17 ± 
0.02 and the expected cross section times efficiency was 0.20 ± 0.02 pb. 

The background included contributions from Z + jet(s) production where one of the jets mimicked an electron or 
photon, multijet production where more than one jet was misidentified as a photon or electron, and ttj production 
followed by decay of each t to ev^Vr. 

Processes where jets mimicked photons, jets mimicked electrons, and double and triple fakes contributed to 
the QCD background. The background and its i?T-dependence were estimated by counting the number of 
ee + jet(s) and 67 + jet(s) events, with the electrons and photons passing the signal cuts and with jet trans- 
verse energy above 10 GeV and 25 GeV, respectively. The probabilities for jets to mimic EM objects were 
determined with a procedure similar to that described in Section V C and observed to be approximately Et 
independent (see Fig. |l^). Table XX contains the probabilities for a jet to mimic photons and electrons 
in the CC and EC, and the final probabilities with the direct photon contribution removed. Multiplying 
these probabilities by the number of jets in these samples led to a background of 0.43 ± 0.06 QCD events. 



Type of Fake 


CC 


EC 


Avg. after direct 7 correction 


Jet^ 7 


(0.84 ±0.08) X 10"^ 


(0.90 ±0.11) X 10"^ 


(0.65 ±0.18) X 10"^ 


Jet— > eTightii 


(0.62 ±0.07) X 10"^ 


(1.5 ±0.2) X 10"^ 


(0.84 ±0.10)) X 10"^ 


Jet > eLoosci 


(1.7±0.1) X 10"^ 


(1.6 ±0.2) X 10"^ 


(1.5 ±0.2) X 10"^ 



TABLE XX. Probability for a jet to mimic a photon or electron, averaged over _Et, for the Z7 ee7 and Z7 fifi-y 
electron and photon selection criteria. 



The Et spectra of the jets allowed the background to be calculated as a function of Et- This is shown in Fig. p2|(a). 
Z boson events where an electron was misreconstructed as a photon and a jet was misreconstructed as the lost 
electron contributed to a bump in the fake photon at Et ^ 50 GeV. Thus, the shape of the fake photon sample was 
parameterized with an exponential function plus a Gaussian. The fit is also shown in Fig. |3^ (a). 

The rr7 background was estimated with a sample of ISAJET events passed through the D0GEANT detector sim- 
ulation and the offline reconstruction algorithm. The total ttj background increased with increasing anomalous 
couplings because more Zj ttj events would have been produced along with the Z — > ecy events. After nor- 
malization with the production cross section and r branching fractions, the expected fraction of the cascading tau 
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decays in the final Z7 sample was /e = (0.10 ± 0.05%), where the uncertainty came from the expected difference 
in acceptance from using Z — > rf ^ eevvvv Monte Carlo to simulate a background which included photons which 
radiated from a charged lepton, as well as the uncertainty in cross section and branching ratios. For the SM couplings 
the rf7 background was negligible. 

To summarize, four — > ee7 candidates were observed. The total background expected was 0.43 ± 0.06 events. 
This corresponds to an observed signal of 3.57l^;9f ± 0.06 events, where the first uncertainty is statistical and the 
second is the uncertainty in the background. The observed signal agrees with the SM prediction of 2.8 ± 0.3 ± 0.2 
events, where the first uncertainty refiects systematics of the Monte Carlo model and the second is the uncertainty in 
the luminosity. 



B. The [ipu channel 

The /Lt/X7 sample was selected from events which satisfied the MUJELE trigger described in Table |. This data 
set corresponded to an integrated luminosity of 13.7 ± O.Tpb"^. At least two muons and one photon were re- 
quired in the event. One muon was required to satisfy the "Tight III" requirements and the other to sat- 
isfy the "Loose 11" requirements of Table ||. It was required that > 15 GeV/c and > 8 GeV/c, 
where \i\ and [12 are the higher and lower muons respectively. The photon, satisfying the "Loose" require- 
ments, was required to have Et > 10 GeV and to be separated from both muons by AT?,^-^ > 0.7, as in the 
electron channel. Two candidates for Zj fj.fj.'y passed these selection criteria. Details on the characteris- 



tics of the candidates and event displays are in Ref. I41f|. Table XXI indicates the number of events which 



survived the last few selection criteria, after the particle identification and kinematic selection were applied. 



Selection Criteria 



No. of Surviving Events 



Particle ID and 
Kinematic Selection 
Trigger Criteria 
A7^u-v > 0.7 



TABLE XXI. Number of fifi^ candidates which passed the selection criteria. 

The efficiencies were calculated as a function of ZZ^ and Zjj couplings with the event generator of Ref. p^ ] 
combined with the parameterized detector simulation. MRSD— ' structure functions [ pTf were used in the event 
generation and the cross section was scaled by a k-factor of 1.34. The level 1 muon trigger efficiency (with two 
chances to trigger on each event), the photon trigger efficiency curve shown in Fig. |3^, and the efficiencies for particle 
identification criteria as discussed in Section ||, were included. The detector acceptance was 20% for SM Z7 /i/i7 
production. The overall efficiency for SM production, for events satisfying the kinematic criteria, was 0.06 ± 0.01. 
The cross section times efficiency for SM production was 0.17 ± 0.03 pb. The efficiency increased with anomalous 
couplings because the muons became more central, increasing their acceptance. 

The background consisted of /i/x -I- jet(s) events, where the jet was misidentified as a photon, and Z — > rr cascade 
decays including a final state photon. The backgrounds where a jet mimicked a photon included Drell-Yan production 
with associated jets, Z -I- jets production and cosmic ray muons (already small because of the tight muon identification 
criteria) in coincidence with jet events. The contribution to the QCD background from bb production was expected 
to be negligible due to the muon isolation requirements and the muon pr threshold. The QCD background was 
estimated from a data sample containing a pair of muons satisfying the same muon identification criteria as the signal 
sample. This sample contained all of the + jet backgrounds in the same proportions as the signal sample. The 
procedure was to count the number of jets of Et > 10 GeV, in events which pass the selection criteria, and multiply 



that number by the probability for a jet to mimic a photon from Table XX . Because different triggers were used in 
collecting the signal and background samples, a scaling factor 1.3 ±0.3 was necessary to account for differences in the 
trigger efficiency and integrated luminosities. Fifteen jets passed the selection criteria; eight jets were in the CC and 
seven in the EC. The result was A^gcn ~ ^■^'^ ^ ^-^^ expected background events. The pr spectrum of the jets in the 
background sample is shown in Fig. p2(b). The fit made to the Et spectrum of the electron channel fakes was used 
to parameterize the Et spectrum of the fakes in the muon channel, with the appropriate normalization, because of 
the much higher statistics of the former. 

The TTj background was estimated using the same procedure as was used for the electron channel. The expected 
fraction of the cascade r decays in the final Z-f ^fx^f sample was 1.4±0.5%. For SM couphngs, the tt^ background 
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was 0.03 ± 0.01 events. 

To summarize, two Z7 fi/ij candidates were observed. The total background expected was 0.05 ± 0.01 events. 
This corresponds to an observed signal of l.QSl^'ag ^ O-^l events, where the first uncertainty is statistical and the 
second is the uncertainty on the background. The observed signal agrees with the SM prediction of 2.3 ± 0.4 ±0.1 
events, where the first uncertainty reflects systematics of the Monte Carlo model and the second is the uncertainty in 
the luminosity. The photon Et spectrum for the combined ee and mi data, expected signal and expected background 
are shown in Fig. 



C. The w channel 



The Z + ^ vv^ signature was a single photon which recoiled against the I^t of the unmeasured neutrino pair. The 
nature of the backgrounds for this channel was very different from the electron and muon channels in that they were 
larger and included contributions from sources to which the previous channels were immune. One background resulted 
from unreconstructed cosmic ray and Tevatron beam related muons which deposited energy in the electromagnetic 
calorimeter through bremsstrahlung. A second important background, occurring at moderately high EtIj), came 
from W ^ ev where the electron was misidentified as a photon due to a missing track. These backgrounds forced 
the analysis to use much stronger particle identification criteria and tighter kinematic selection than the analyses 
presented above. 

The candidate sample was selected from events passing the ELE-HIGH trigger of Table |. A selection on the "event 
quality" removed events with noisy cells in the calorimeter, second EM objects with Et > 5 GeV, or when the 
calorimeter was recovering due to a large pulse from Main Ring associated energy deposition. Both the photon Et 
and the I^t were required to be greater than 40 GeV to reduce the background from W boson decays. Events with 
muons in the central region were rejected to reduce cosmic ray backgrounds. Events with jets of Et greater than 
15 GeV were also rejected; by limiting the pt boost of the events, the kinematic range of I^t and E^ from the W 
boson background was reduced. The strictest photon requirement, "Tight" , was used to reduce the backgrounds 
from cosmic rays, from beam related muons, and from W boson decays. Lastly, the calorimeter was searched, in a 
road about the line defined by the vertex located by the EMVTX package (see Section [I C 9 ) and the energy- weighted 
center of the photon shower in the CC, for energy deposition consistent with the passage of an unreconstructed cosmic 
ray which might have radiated the photon. This algorithm (MTC) tracked the muon energy deposition through the 
longitudinally and azimuthally segmented towers of the calorimeter. 

Four events remained after all selection criteria were applied. The photons in the four events had Et of 
41, 41, 46, and 68 GeV. Table XXII shows the number of events remaining after each of the selection criteria. 



Selection Criteria 


Number of Events Remaining 


Trigger, Event Quality and Kinematics 


1887 


77^ < 1.0 or 1.5 < \rij\ < 2.5 


1637 


"Loose" Photon Criteria 


1448 


<^ = 


1098 


No Jet with Et > 15 GeV 


480 


"Tight" Photon Criteria 


5 


MTC selection 


4 



TABLE XXII. Event selection criteria for the Z"/ wy analysis. 



The ELE-HIGH trigger was completely efficient for Et{j) > 40 GeV. The efficiency for the photon identification 
criteria (excluding the fiducial requirements), the event quality, the muon and jet vetoes, and the MTC selection 
criteria were estimated using the Z ^ ee candidates collected using the same trigger as the signal. The efficiency 
for the jet veto was cross checked with a Next-to-Leading logarithm Monte Carlo generator ^] . The calculated 
efficiency loss agreed with the measurement. Table XXIII contains a summary of these efficiencies. The efficiency for 
the fiducial selection came from the event generator of Ref. pq] combined with the parameterized detector simulation. 
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Selection Criteria 


CC 


Combined 


EC 


Event Quality 




0.981 ± 0.002 




Photon ID Criteria 


0.61 ±0.02 




0.66 ±0.05 


N^^ = 




0.988 ± 0.002 




No Jet with Et > 15 GeV 




0.84 ± 0.02 




MTC Selection 


0.97 ±0.02 






Total 


0.48 ± 0.02 




0.54 ±0.04 



TABLE XXIII. Efficiency, excluding photon fiducial requirements, for the analysis. The uncertainties indicated 

are statistical only. 



The cosmic ray and beam halo backgrounds, due to unreconstructed muons which radiated a photon as they passed 
through the calorimeter, was estimated using cosmic ray and beam halo muons identified in the data. 

Two samples of this kind of background event which radiated a photon into the CC were identified. The first 
sample was identified by applying all the selection criteria except for the HITSINFO criteria of the "Tight" photon ID 
requirements, the CF muon veto, and the MTC requirement. The event was required to have a reconstructed muon. 
The rejections for the HITSINFO and the MTC selection criteria were determined from this sample. The reconstruction 
efficiency for the muon background (ecosmic ^l) were estimated from a second sample of events, dominated by cosmic 
ray muons, which passed the same selection criteria as the former sample (excluding the requirement that the muon 
was reconstructed) and failed the EMVTX criteria. The inefficiency was 1 — Ccosmic ^ — 0.34 ± 0.03. The background 
to Z"f was then determined from the number of events in the former sample, modified by the rejection provided by 
the HITSINFO and MTC criteria and by the factor (1 — £cosmic fj.) / ^cosmic ti- The resulting expected background was 
1.4 ± 0.6 events for CC photons. 

For the muon background events with a photon in the EC, a two-sample procedure analogous to that described 
above for the CC was used to estimate the background due to beam-halo muons. It was measured to be 0.38 ± 0.23 
events, where the large fractional uncertainty is due to the low statistics of one of the tagged samples. 

The muon background rejection was found to be independent of the photon Et for both the CC and EC regions, 
so the background events were used to make a parameterization of the photon spectrum. The doubly-identified 
background muons were used for the CC. For the EC, the larger of the two available background samples was used. 
The result of the fit is shown in Fig. ^ together with the background data. 

Another background comes from W ^ ev events where the electron is misidentified as a photon due to tracking 
inefficiency. The kinematic requirements combined with the jet veto rejected some of the misidentified W boson 
decays, but additional rejection was required. It came from the HITSINFO criteria of the "Tight" photon selection. A 
procedure similar to that used to estimate the background due to muon bremsstrahlung was used in this case. A nearly 
pure sample of tagged W ^ ev events with and Ij^T both greater than 25 GeV was obtained by applying all of the 
event selection criteria except for the HITSINFO and MTC (CC only) criteria and by requiring an electron with a good 
track match significance, TMS < 10 (see Section II C 5| ) , instead of a photon. The rejection, Rh, due to HITSINFO 
came from a sample oi W ev events which failed the electron tracking requirement. Including small corrections to 
account for the fraction of these mistracked W events, lost from this sample because of an overlapping random track, 
and to account for cosmic ray bremsstrahlung in this background sample, Rh was found to be 48 ± 12 (43 ± 14) for 
the CC (EC). With the measured efficiency of the track finding, the track match significance and the MTC criteria 
(ct, Cct, and eMTC, respectively), the W ^ ev background was simply the number of tagged W ^ ev events times 
(1 - €T)/{eT eMTC Rh)- This was found to be 2.2 ± 0.6 (1.8 ± 0.6) events for the CC (EC). The Et spectrum of 
the expected background is shown in Fig. ^ Parameterizations for the expected E^ spectrum of background were 
derived from the W ^ ev events. They are shown in Fig. 

The possible QCD backgrounds included: multijet production, where a jet was misidentified as a photon and the 
l^T results from mismeasurement of a jet or from neutrinos in a jet; direct photon production (jet + photon) where 
a jet contributes to I^t'i and Z + jets — > vD + jets where a jet is misidentified as a photon. However, the size of the 
QCD backgrounds fell rapidly as the photon Et and thresholds were raised. The backgrounds were found to be 
negligible for E^ and > 35 GeV. 



The total background was 3.6 ± 0.8 (2.2 ± 0.6) for the CC (EC). It is summarized in Table |XXIV 
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Background 


atCC 
'-^hck 




Muon Bremsstrahlung 


1.4 ±0.6 


0.38 ±0.23 


W -^eu 


2.2 ±0.6 


1.8 ±0.6 


QCD Sources 


negligible 


negligible 


Total 


3.6 ±0.8 


2.2 ±0.6 



TABLE XXIV. Number of expected background events in the CC and EC for the Z'y wy analysis. 



The expected numbers of signal events for SM and for anomalous couplings were estimated using the leading order 
event generator of Ref. jl^ combined with the parameterized detector simulation, including a pt{Z) spectrum from 
the Z boson data to mimic the effects of the jet veto on the acceptance. The energy scale for the underlying event 
was determined by comparing the pt as determined from the electrons and hadronic recoil in low px Z ee events. 
As in the charged lepton analyses, the cross section was scaled by a fc-factor of 1.34 and the MRSD— ' structure 
functions [ pTj were used in the event generation. A 12% uncertainty in the cross section resulted primarily from 
the choice of parton distribution functions, modeling of the jet veto, modeling of the detector, and the detector 
efficiency. Table XXV presents a summary of the expected signal and background as well as the number of events 
seen with photons in the CC and EC. The SM signal was expected to be 1.8 ± 0.2 events with a 5.8 ± 1.0 event 
background. For comparison, the expected number of signal events for anomalous couplings was approximately a 
factor of 9 higher for h^Q = 3, /ifg = 1. Four candidates were observed, consistent with the SM expectations. The 
photon spectra expected for the signal and background, as well as that seen in the data are shown in Fig. |3^. 



Region 


NsM 


Nbck 


Data 


CC < 1.1) 


1.4 ±0.2 


3.6 ±0.8 


3 


EC (1.5 < |77^| < 2.5) 


0.39 ±0.05 


2.2 ±0.6 


1 


Total 


1.8 ±0.2 


5.8 ± 1.0 


4 



TABLE XXV. Number of expected uvy events assuming SM couplings, number of expected background events, and observed 
signal. 



D. Limits on Anomalous ZZ"/ and ^77 Couplings 



To set limits on the anomalous coupling parameters, the observed Et spectrum of the photons in the three 
channels was fit with the MC predictions plus the estimated background (summarized in Tabic XXVI). The 
binned likelihood method described in Appendix 2 was used. To exploit the prediction that anomalous cou- 
plings lead to an excess of events with high Et photons, a high Et bin with no events was used in the 
fit. Common systematic uncertainties, including photon identification efficiency, integrated luminosity, choice 
of parton distribution fimctions, and choice of pt{Z) distribution were treated as discussed in Appendix 2. 



Channel 


ee 




uv 




2.8 ±0.3 ±0.2 


2.3 ±0.4 ±0.1 


1.8 ±0.2 




0.43 ± 0.06 


0.05 ±0.01 


5.8 ± 1.0 


j\^observed 


4 


2 


4 



TABLE XXVI. Summary the of Z^ — > dileptons analyses including the number of SM events expected, expected backgrounds, 
and number of candidates observed. 



The form-factor scale dependence of the result was studied. The chosen value of A = 500 GeV was close to the 
sensitivity limit of the experiment for /120 and /140 for the ee-\-\x\x channels; for larger values of A partial wave unitarity 
was violated for certain values of the coupling parameters allowed at 95% CL. With the vv and combined analysis, A 
could be extended to 750 GeV without violating unitarity. In that case, tighter limits on anomalous couplings could 
be obtained. 

Figures pq to show the coupling limits for the CP-conserving ZZ7 parameters. The shapes of the ^'77 
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limit contours were similar. 
CL. Figure p 



Figure |3^ shows the results of the fit for the ee and ^/z channels at 68% and 95% 
shows the results of the fit for the vv channels and for the three channels combined at 95% CL. 
The form factor scale A — 500 GeV was used in these two figures. Figure 1^ shows the 95 % CL h mits for 
A = 750 GeV. The 95% confidence level limits on CP-conserving couplings are given in Table XXVII, Shown 
are the limits for the ee + p^ii channels, the vv channel, and the limits from the three analyses combined. Lim- 
its on the CP- violating couplings were numerically identical to the corresponding CP-conserving couplings with 
the single exception that —0.86 < /ifg < 0.87 (to be compared with —0.87 < ft-fg < 0.87) for the vv^ analysis. 



ee + /i/i 


vv 


Combined Limits 


Limits with A = 500 GeV 


-L8 < hia < 1.8 


-0.87 < hia < 0.87 


-0.78 < hio < 0.78 


-L9 < < 1.9 


-0.90 < /iJo < 0.90 


-0.81 < ft]o < 0.81 


-0.5 < h^o < 0.5 


-0.21 < hio < 0.21 


-0.19 < hfo < 0.19 


-0.5 < /i]o < 0.5 


-0.22 < /i]o < 0.22 


-0.20 < hJo < 0.20 


Limits with A = 750 GeV 




-0.49 < /ifo < 0.49 


-0.44 < fefo < 0.44 




-0.50 < /i]o < 0.50 


-0.45 < hig < 0.45 




-0.07 < hfo < 0.07 


-0.06 < h^o < 0-06 




-0.07 < hjg < 0.07 


-0.06 < hjg < 0.06 



TABLE XXVII. Limits on CP-conserving ZZy and ^77 anomalous coupling parameters for the ee + fijj,, vv and combined 
analyses. These axes hmits are at 95% confidence level with A = 500 and 750 GeV. 



X. CONCLUSIONS 



Four gauge boson pair production processes and corresponding trilinear gauge boson coupling parameters were 
studied using the data from 1.8 TeV pp collisions collected with the D0 detector during 1992-1993 Tevatron collider 
run at Fermilab. The data sample corresponded to an integrated luminosity of approximately 14 pb~^. 

Searches were made for deviations from the SM. This would have been manifest as an enhancement in the production 
cross section and Et spectrum of the bosons. In the analyses of the VL7 final states, the photon Et spectrum 
was compared with the expectations of the SM and used to produce limits on anomalous WW^ couplings. A 
limit on the cross section for WW — > dilepton led to a limit on anomalous WW-f and WWZ couplings. In the 
WW/WZ e:^ jet jet analysis, the Et spectrum of the W and Z bosons was used to produce limits on anomalous 
WW'y and WWZ couplings. The 11^7, WW and WW/WZ analyses were combined to produce limits on WW^ 
and WWZ couplings. Finally, in the analysis of Z7 final states, the photon Et spectrum was compared to the 
expectations of the SM and used to produce limits on anomalous ZZ7 and Z77 couplings. No deviations from the 
SM were observed. 

The W^ analysis yielded 23 candidate events where the W boson was identified by its leptonic decay products, a 
high "pt electron (11 events) or muon (12 events), and a neutrino inferred by large I/t in the event. The expected 
backgrounds for the electron and muon channels were 2.0±0.9 and 4.4±1.1 events, respectively. Using the acceptance 
for the SM W^ production events, the W^ cross section (for photons with Elp > 10 GeV and ATZi^ > 0.7) was 
calculated from the combined e + fi sample to be a{Wj) = ISS't^sistat) ± 21{syst) pb. A binned maximum likelihood 
fit was performed on the Pj. spectrum for each of the W{ev)j and W{fj,v)j samples to set limits on the anomalous 
coupling parameters. The limits on the CP-conserving anomalous coupling parameters at the 95% CL were —1.6 < 
An^ < 1.8 (A^ = 0), -0.6 < < 0.6 {Ak^ = 0) using a form factor scale of A = 1.5 TeV. The U{l)EM-only 
couphngs of the W boson to a photon, which lead to k.^ = and A.^ = 0, and thereby fiw — e/2mw and = 0, 
was excluded at the 80% CL, while zero magnetic moment {^iw — 0) was excluded at the 95% CL. 

The search for WW events where both of the W bosons decay leptonically to ev or iiv yielded one candidate 
event with an expected background of 0.56 ± 0.13 events. The upper limit on the cross section for SM W boson 
pair production was estimated to be 87 pb at the 95% CL. The limit on the cross section was translated into 
limits on the anomalous coupling parameters. The limits on the CP-conserving anomalous coupling parameters were 
-2.6 < Ak < 2.8 (A = 0) and -2.1 < A < 2.1 (Ak = 0) at the 95% CL, using a form factor scale of A = 900 GeV, 
where Ak = Akj = Akz and A = A-y — A^ were assumed. 
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The analysis of the WW and WZ production events in the electron + jets channels, where one W boson decayed 
into ev and the second W boson or Z boson decayed into two jets, yielded 84 candidate events with an expected 
background of 75.8 ± 13.3 events, while the SM predicted 3.2 ± 0.6 signal events. A maximum likelihood fit was 
performed on the px spectrum of the W boson, computed from the Et of electron and the I^t, to set limits on the 
anomalous couplings. The limits on the CP-conserving anomalous coupling parameters were —0.9 < Ak < 1.1 (A = 0) 
and —0.6 < A < 0.7 (Ak = 0) at the 95% CL, using a form factor scale of A = 1.5 TeV, where Ak = Ak-^ = Ak^ and 
A = A-y = A^ were assumed. 

The W^ production process is sensitive only to the WW'y coupling parameters. WZ production is sensitive only to 
the WWZ couplings. On the other hand, W pair production is sensitive both to the WW^ and the WWZ coupling 
parameters. Using assumptions on the relationship between the WW'y and WWZ coupling parameters, these three 
analyses were combined to set the tightest limits on the coupling parameters. A maximum likelihood fit was performed 
on the three sets of data simultaneously using a common form factor scale of A = 1.5 TeV and the assumption that 
Ak = Ak^, = Ak^ and A = A-y = Xz- The limits obtained at the 95% CL were -0.71 < Ak < 0.89 (A = 0) and 
—0.44 < A < 0.44 (Ak = 0). These are the tightest limits on WW-f and WWZ couplings presently available. 

Zj final states in pp collisions are produced from the Z boson-quark and photon-quark couplings in the SM. The 
ZZj and Zjj couplings, which can produce Zj final states, are absent in the SM. The Z7 analysis yielded a total of 10 
candidate events; 4 events with Z — > ee, 2 events with Z —^ f^ii and 4 events with Z —^ vv. The expected backgrounds 
for the ee and ///i channels were 0.43 ± 0.06 events and 0.05 ± 0.01 events, respectively. The expected background for 
the vv channel was 5.8 ± 1.0 events. The sum of the SM prediction and the expected background was 2.8 ± 0.3 ± 0.2 
(2.3 ± 0.4 ± 0.1) events for the electron (muon) decay modes and 7.6 ± 1.0 events for the neutrino decay mode. A 
maximum likelihood fit was performed on the pt spectrum of the photons to set limits on the anomalous coupling 
parameters. The 95% CL axes limits on the CP-conserving coupling parameters are —0.78 < h^^ < 0.78 (/i^q — 
0);-0.19 < /ifo < 0.19 (/ifo = 0),-0.81 < hj^ < 0.81 {hj„ = 0);-0.20 < hj^ < 0.20 {hj^ = 0), using a form factor 
scale of A = 500 GeV and -0.44 < /ifg < 0.44 (/ifg = 0); -0.06 < hf^ < 0.06 (/iff, = 0), -0.45 < hj„ < 0.45 {hj^ = 
0); —0.06 < hjg < 0.06 (^130 = 0), using a form factor scale of A = 750 GeV. The limits obtained in this measurement 
are the most stringent limits on anomalous ZVj couplings currently available. 

All of these limits on the anomalous coupling parameters will be significantly improved when the analyses of the 
data taken during the 1994-1995 Tevatron collider run, which corresponded to approximately 80 pb~^, are completed. 
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XI. APPENDIX 1: TRILINEAR GAUGE BOSON COUPLING PARAMETERS 



A. WW'y and WWZ coupling parameters 

The tree-level Feynman diagrams for qq — > Wj, qq — > WW and qq — s- WZ production processes are shown in 
Figs. |l|, H and ^. A formalism has been developed to describe the WW'y and WWZ vertices for the most general 
gauge boson self-interactions [f 15|. The Lorentz invariant effective Lagrangian for the gauge boson self-interactions 



contains fourteen dimensionless coupling parameters, seven each for WW^ and WWZ: 

Cwwv/gwwv = igX [wl^w^^V' - WjKW^^") 
+iKvWlW,V^'^ + ^^wlW|:V''^ 
-g^WlW^id^V + d''V'') 

+ikvwlw,w^ + '-^wlw;^v''\ 
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where VF'' denotes the field, W^, = d^W. - d,W^, V^, = - 9.^^, %u = ^e^-P"^^", and {A B) = 
A{d^B) — {dfj,A)B, V = J, Z and is the mass of the W boson. The overaU couphngs gwwv are gwwj = — e and 
gwwz = — e(cot0u,) as in the SM, where e and 9^ are the positron charge and the weak mixing angle. The couplings 
Xv and Ky conserve C and P. The couplings gY are odd under CP and C, are odd under C and P, and Ky and 
Ay are odd under CP and P. 

In the SM, all the couplings are zero with the exception of gY and Ky {gj = gf = k-^ = = 1). Electromagnetic 
gauge invariance restricts gj, gj and gj to the SM values of 1, 0, and 0. The SU{2) l x t/(l)Y gauge invariance requires 
A = A-y = Az and A = A-^ = Az- If the photon and Z boson couplings are assumed to be equal, then = ffi = 1 and 
5f = ff4 = fff = ff5 = 0. 

With the non-SM coupling parameters, the amplitudes for the gauge boson pair production grow with energy, 
eventually violating tree-level unitarity. Using dipole form factors for anomalous couplings, 

Ak(s) = 



(1 + J/A2)2 

with a form factor scale. A, the unitarity is restored. The scale. A, is constrained by 

6.88 



A< 



(k-1)2 + 2A2 + 2A2 



-,1/4 

TeV, (11.1) 



if the photon and Z boson couplings are assumed to be equal. 

The CP-conserving WW^ coupling parameters are related ||] to the magnetic dipole moment (iiw) and electric 
quadrupole moment (Q^) of the W boson: 



The CP-violating WWj coupling parameters are related to the electric dipole moment (dw) and magnetic quadrupole 
moment (Q^y) of the W boson: 



The CP-violating WWj couplings A-,, and are tightly constrained by measurements of the neutron electric dipole 
moment to \kJ, \\J < 10^^ H. 



B. ZZ^ and ^77 Couplings 

Theoretical calculations of the tree-level cross section for Z7 production for SM and anomalous couplings have been 
performed [Q. Assuming only electromagnetic gauge invariance and Lorentz invariance, the vertex function for the 
ZZj and Zjj interaction can be described with the following form: 

CvihYiq^g-P-qf) 
+ ^P'^{{p.q^)g^^f^-q^p0) 

+ hYe^'"'''q2p + ^P''ef'^f"'Ppq2a), 

where V indicates a photon or Z boson, Cz is (P^ — (^IM\ and C^ is P"^ /M"^^ qi and 52 are the momenta of the 
outgoing particles and P is the momentum of the virtual boson. An overall normalization factor oi gzz'^ = gz-yy = e. 
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left out in the equation, is used. These couplings are C odd dimensionless functions of gj, ^'^'^ ^ 5 i-^- s. In 
addition, hX and are P even, and thus violate CP. The other pair, and /i^^, are CP conserving. 

In order to avoid violating ^-matrix unitarity, the couplings should asymptotically approach zero (their SM value) 
at high energies Therefore, the ZV^ couplings have to be energy-dependent and are thus modified with form 

factors hY {q\, q^, P^which vanish at high g^, g|, or P^. However, since q\ ~ m^, « and P^ = s, only the high s 
behavior should be included in the form factor for the qq — > diagrams. The convention is to use a generalized 
dipole form factor such that: 



hY{ml,Q,s) 



hi 



(1 + S/A2)- 



The constraints on the HYq can be derived from partial wave unitarity of the general // process |3l| , [45 

Assuming only one coupling is non-zero at a time, the following unitarity limits can be derived for A ^ mz |1§[| , [46 



(|n)" 0.126/0.151TeV^ 



10 hl-30 I ^ (2„„i)„-3/2 A3 



z/,. .^z/-r. ^ (f^)" 2.1/2.5 -lO-^TeV^ 



From the above equations, unitarity is satisfied for n > 3/2 for hY^, and n > 5/2 for hY^. In this paper n = 3 for hY^, 
and n = 4 for hY 4 are used. This choice ensures the same asymptotic energy behavior for the hY^ and hY 4 couplings. 
The dependence of results on the choice of n is discussed in [Q . 

The anomalous couplings hY are related to the Zj transition dipole and quadrupole moments. The CP-even 
combinations of hY and correspond to the electric dipole and magnetic quadrupole transition moments; the CP- 
odd combinations of hY and hY correspond to magnetic dipole and electric quadrupole transition moments. The 
relations between the couplings and moments depends on both the center of mass energy a/I and on the momentum 
of the final state photon [^. They are: 

J ^ ^ f z 1 z \ 

Q% - ^VlO(2/ifo), 

Q", = ^\/l0(2/ifo), 

where dz^ ifJ-Zr) is the transition electric (magnetic) dipole moment, Q^t iQ"zT^ the transition electric (magnetic) 
quadrupole moment, and k is the photon energy. 

XII. APPENDIX 2: BINNED LIKELIHOOD FIT 

A binned likelihood fit was applied to the px spectra of 7 and W to set limits on the anomalous coupling parameters. 
The observed numbers of events (iV^) in a particular pT bin can be described in terms of the numbers of expected 
signal events {rii) and background events (6^) using a Poisson distribution function: 

where and in some cases hi are functions of the anomalous coupling parameters. 

The uncertainties in the rii and hi were incorporated by convoluting with Gaussian distributions: 
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p' 



oo oo 



-oo — oo 



^ 2,7? 



27r (TnCTb 



where and /& are multiplicative factors to rii and bi with mean values of 1.0; (t„ and tJb are the fractional uncertainties 
of Hi and bi . These uncertainties include the uncertainties in the integrated luminosity and the theoretical prediction 
of the signal and background cross sections. To exploit the prediction that anomalous couplings lead to an excess of 
events with high Et photons or jets (depending on the analysis), a bin with no events at high Et was used in the 
fit. The bin boundary was selected sufficiently above the highest observed transverse momentum event in the data 
sample that the detector resolution could not move the last data point across the boundary. For more detail, see [^ . 
The joint probability of all px bins is then 



n ^/ 



i=i 

where Nmax is the number of pr bins. The log likelihood function of this joint probability is defined as 

C = -ln(P). 

The limits on the coupling parameters were obtained by maximizing this quantity. The 95% confidence level limit 
on the parameters of the log-likelihood function (the coupling parameters for the case here) is the contour where the 
log-likelihood is 1.92 lower than the maximum. 
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FIG. 1. Leading-order Feynman diagrams for vector boson pair production. The assignment of Vb, Vi, and V2 depends on 
the final state: Wj, WW, WZ or Zy. 
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FIG. 2. The leading-order radiative Feynman diagram for production where the photon is the result of bremsstrahlung 
from a final state lepton. 
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FIG. 3. The predicted Et spectrum of photons in the production for SM and anomalous WW^ coupUngs. Radiative 
diagrams are included. The requirement that AR^-y > 0.7 has been made; otherwise the cross section diverges at low E?^. 
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FIG. 4. The predicted cross section for WW production as a function of anomalous coupling parameters A and A/t, assuming 
the WWZ and WWj couplings are equal, with A = 1000 GeV. 
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FIG. 5. The pr spectrum of W bosons in WW/WZ production corresponding to approximately 10 pb ^ of collisions. The 
theoretical assumptions for the anomalous coupling spectra are that A-y = Xz and Ak^ = Akz, with A = 1000 GeV. 
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FIG. 6. The leading-order radiative Feynman diagrams for production where the photon is the result of bremsstrahlung 
from a final state lepton. This decay mode only applies to final states involving charged leptons. 
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FIG. 7. The spectrum, da/dp^, for SM events from final state radiation (dashed line), initial state radiation (dotted 
line), the combination of initial state and final state radiation (solid line) as well as for two anomalous ZZ-y couplings. 
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FIG. 8. Perspective view of the D0 detector. Also shown are the movable support platform, the Tevatron beampipe centered 
within the detector and the Main Ring beampipe which penetrates the muon system and calorimeter above the detector center. 
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FIG. 9. The D0 central tracking detector system. 
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FIG. 10. The D0 calorimeters. 
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FIG. 11. Muon selection variables including (a) the bend view impact parameter, (b) the nonbend view impact parameter, 
(c) the energy in the calorimeter tower plus the nearest neighboring towers around the muon, and (d) the to resulting from the 
track fit with the muon time-of-origin as a parameter. 
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FIG. 12. Electron selection efficiencies for (a) the electron selection efficiency as a function of the fraction of the energy 
deposited in the EM calorimeter, (b) the H-matrix x^, (c) the isolation variable fiso, and (d) the track-match-significance, 
TMS. The solid circles are for CC electrons and the open diamonds are for EC electrons. 
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FIG. 13. The efficiency of the "Loose" photon H-matrix selection criteria as a function of 
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FIG. 14. (a) the dE/dx measured in the CDC for electrons from Z boson decays, (b) the dE/dx measured in the CDC for 
EM clusters in an inclusive jet sample. 
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FIG. 15. (a) jy-dependent probability for photon conversion (Pc) in the material in front of the CDC and FDC ; (b) 
77-dependent efficiency of the photon identification (tj) for high pr photons for the Z7 "Loose" selection criteria. The un- 
certainty shown includes the statistical uncertainty plus a common systematic uncertainty of 5%. 
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FIG. 16. Probability distributions Pxy and Prz- (a) and (b) are from photons resulting from cosmic ray bremsstrahlung. 
(c) and (d) are from electrons from Z boson decays. 
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FIG. 17. The probability of a jet to be misidentified as a photon as a function of Et (before the removal of the contribution 
from direct photons) for the "Loose" photon selection criteria in the CC (a) and EC (b). 
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FIG. 18. Distribution of (a) E^, (b) ARe^ and (c) Mri-fi;!^) for the W{eu)-y + W{fj,i')j combined sample. The points 
are data. The shaded areas represent the estimated background, and the sohd histograms are the expected signal from the 
Standard Model plus the estimated background. 



55 




FIG. 19. Limits on (a) CP-conserving anomalous coupling parameters Ak and A, and on (b) the magnetic dipole, nw, and 

electric quadrupolc, Qiy, moments. The ellipses represent the 68% and 95% CL exclusion contours. The symbol, •, represents 
the Standard Model values, while the symbol, indicates the l7(l)£M-only coupling of the W boson to a photon, Ak = — 1 
and A = {nw = e/2mw and Ql^r = 0). 
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FIG. 20. AcI>{pt,1?t) vs I^t distributions for WW {ev){ev) with the SM couplings, Z ^ ee, Z ^ tt ^ eevvvv and 
WW — > eevv) with the non-SM coupUngs. 
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FIG. 21. Ei^'^'^ for Monte Carlo WW (open histogram) and tt events (shaded histogram) with M^p = 160 GeV/c^ (/ Ldt ~ 20 
fb"^). Events with B^*^ > 40 GeV were rejected. 
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FIG. 22. A<j>(p^, I^t) vs distributions for WW — > e/i with the SM coupUngs, Z — »• rr — » e^.vi'i'i' and tyW^ eji with 
non-SM couplings. 
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FIG. 23. .^J vs Acj>{pf^\^T) distributions for WW — » nni/v with SM couplings, Z — » nfj,, and tt events. 
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FIG. 25. The distribution of the transverse mass of the electron and for the data (points), major backgrounds (soUd 
Une), 10 times SM WW signal (dotted), and WW with Ak = 2 A = (dot-dash). The backgrounds are normalized as described 
in the text. 
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FIG. 26. The distribution of the dijet invariant mass for the data and major backgrounds. The backgrounds are normalized 
as described in the text. The arrows indicate the region accepted by the dijet mass selection criterion. 
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FIG. 27. Efficiency for reconstructing the dijets and for the dijet mass selection for W —^ jj vs. pt(W). The solid crosses 
are the results from ISAJET. The dashed crosses are results from PYTHIA. 
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FIG. 28. Distribution of the for the WW/WZ candidates (sohd) and the QCD fake sample (dashed) before the dijet 
mass selection. 
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FIG. 29. pt distributions of the ef system: data (points), W+ > 2 jets background (dotted), total background (solid), and 
Monte Carlo predictions for the SM (dashed) and non-SM couplings Ak = 2 A = (dot-dashed) WW production. 
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FIG. 30. Contour limits on anomalous coupling parameters at the 95% CL (inner curves) and limits from S'-matrix unitarity 

(outer curves) for the assumptions a) A-^ — Xz and Anj = Akz, b) HISZ relations, c) SM WW^ couplings and d) SM WWZ 
couplings. A = 1500 GeV was used for (a), (b), and (c). A = 1000 GeV was used for (d). 



67 




FIG. 31. 95% CL limits (inner contour) on A and A/t, assuming the WW'y and WW Z couplings are equal and A = 1500 
GeV, from the combined W^, WW, and WZ results. The outer contour is the limit from s-matrix unitarity. 
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FIG. 32. a) QCD background in the ee7 channel. The bump around 50 GeV/c is due to the electrons from Z boson decays 
where a jet mimics an electron and an electron mimics a photon, b) QCD background in the /^/i7 channel. The fit is described 
in the text. 
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FIG. 33. Efficiency as a function of photon Et for the level 1 EM trigger with threshold at 7 GeV. The minimum allowed 
photon Et is 10 GeV. 
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FIG. 34. Transverse energy spectrum of photons in 667 and fij^f events. The shadowed bars correspond to the data, the 
hatched curve represents the total for background, and the soUd hne shows the sum of SM predictions and background. 
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FIG. 35. Shape of the bremsstrahlung photon spectrum for the background due to cosmic ray and beam halo muons in the 
a) CC and b) EC calorimeters. The solid lines are the resulting fitted parameterizations. 
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FIG. 36. Shape of the photon spectrum for the background due to W ev in the Z'y — > z/1/7 analysis, (a) is the spectrum 
in the CC and EC. (b) and (c) are the individual spectra in the CC and EC where the fits to the background are shown (lines). 
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FIG. 37. The photon Et spectra for the Z"f i^iyy data (points), the background (sohd Une), the expected signal (shaded), 
and the sum of the expected signal and background with uncertainties (dotted). 
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FIG. 38. Limits on the CP-conserving anomalous ZZ'y coupling parameters /ifo and /ifo. The solid ellipses represent 68% 
and 95% CL exclusion contours for the ee and jifx combined analysis. The dashed curve shows limits from partial wave unitarity 
for A = 500 GeV. 
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FIG. 39. Limits on the CP-conserving anomalous ZZ^ coupling parameters Tifo and ftfo- The solid ellipses represent the 
95% CL exclusion contours for the vv and for combined ee, /U/U, and vV analyses. The dashed curve shows limits from partial 
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FIG. 40. Limits on the CP-conserving anomalous ZZ^ coupling parameters /ifo and ftfo- The solid ellipses represent the 
95% CL exclusion contours for the vv and for the combined ee, /U/U, and vv analysis. The dashed curve shows limits from 
partial wave unitarity for A = 750 GeV. 
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